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ABSTRACT 

A Chandra ACIS S3 observation of the X-ray faint elliptical galaxy NGC 4697 resolves much of the 
X-ray emission (61% of the counts from within one effective radius) into 90 point sources, of which 
~80 are low mass X-ray binaries (LMXBs) associated with this galaxy. The dominance of LMXBs 
indicates that X-ray faint early-type galaxies have lost much of their interstellar gas. On the other 
hand, a modest portion of the X-ray emission from NGC 4697 is due to hot gas. Of the unresolved 
emission, it is likely that about half is from fainter unresolved LMXBs, while the other half (^23% of 
the total count rate) is from interstellar gas. The X-ray emitting gas in NGC 4697 has a rather low 
temperature (kT = 0.29 keV). The emission from the gas is very extended, with a much flatter surface 
brightness profile than the optical light, and has an irregular, L-shaped morphology. The physical state 
of the hot gas is uncertain; the X-ray luminosity and extended surface brightness are inconsistent with a 
global supersonic wind, a partial wind, or a global cooling inflow. The gas may be undergoing subsonic 
inflation, rotationally induced outflow, or ram pressure stripping. X-ray spectra of the resolved sources 
and diffuse emission show that the soft X-ray spectral component, found in this and other X-ray faint 
ellipticals with ROSAT, is due to interstellar gas. The cumulative LMXB spectrum is well-fit by thermal 
bremsstrahlung at kT = 8.1 keV, without a significant soft component. 

NGC 4697 has a central X-ray source with a luminosity of Lx = 8 x 10 38 ergs s _1 , which may be due 
to an AGN and/or one or more LMXBs. At most, the massive black hole at the center of this galaxy is 
radiating at a very small fraction (< 4 x 10~ 8 ) of its Eddington luminosity. 

Three of the resolved sources in NGC 4697 are supersoft sources. In the outer regions of NGC 4697, 
seven of the LMXBs (about 20%) are coincident with candidate globular clusters, which indicates that 
globulars have a high probability of containing X-ray binaries compared to the normal stellar population. 
We suggest that all of the LMXBs may have been formed in globulars. The X-ray-to-optical luminosity 
ratio for the LMXBs in NGC 4697 is L X (LMXB, 0.3-10 keV)/L B = 8.1 x 10 i9 ergs s" 1 L^, which 
is about 35% higher than the value for the bulge of M31. Other comparisons suggest that there are 
significant variations (factor of >2) in the LMXB X-ray-to-optical ratios of early-type galaxies and spiral 
bulges. The X-ray luminosity function of NGC 4697 is also flatter than that found for the bulge of M31. 
The X-ray luminosities (0.3-10 keV) of the resolved LMXBs range from ^5 x 10 37 to ^2.5 x 10 39 ergs s _1 . 
The luminosity function of the LMXBs has a "knee" at 3.2 x 10 38 ergs s _1 , which is approximately the 
Eddington luminosity of a 1.4 M & neutron star (NS). This knee appears to be a characteristic feature of 
the LMXB population of early- type galaxies, and we argue that it separates black hole and NS binaries. 
This characteristic luminosity could be used as a distance estimator. If they are Eddington limited, the 
brightest LMXBs contain fairly massive accreting black holes. The presence of this large population 
of NS and massive BH stellar remnants in this elliptical galaxy shows that it (or its progenitors) once 
contained a large population of massive main sequence stars. 

Subject headings: binaries: close — galaxies: elliptical and lenticular — galaxies: ISM — X-rays: 
galaxies — X-rays: ISM — X-rays: stars 



1. INTRODUCTION 

X-ray observations, starting with the Einstein X-ray 
Observatory, have shown that elliptical and SO galaxies are 
luminous sources of X-ray emission (e.g., Forman, Jones, 
& Tucker 1985). At least for the X-ray luminous early- 
type galaxies (defined as those having a relatively high 
ratio of X-ray to optical luminosity Lx/Lb), it is clear 
that the bulk of the X-ray luminosity is from hot interstel- 
lar hot (~10 7 K) interstellar gas (e.g., Forman et al. 1985; 
Trinchieri, Fabbiano, & Canizares 1986). Presumably, the 
material comprising the interstellar medium (ISM) orig- 



inated from stellar mass loss in the galaxy. There is a 
strong correlation between the X-ray and blue luminosi- 
ties of early-type galaxies Lx oc L^ 7-3 (Canizares, Fab- 
biano, & Trinchieri 1987; White & Davis 1997; Brown & 
Bregman 1998). However, there is a large dispersion in the 
X-ray luminosities of early-type galaxies of a given optical 
luminosity. Two galaxies with similar blue luminosities 
might have X-ray luminosities that differ by as much as 
a factor of 100 (Canizares et al. 1987; Fabbiano, Kim, & 
Trinchieri 1992; Brown & Bregman 1998). We will refer 
to galaxies which have a very low Lx /Lb ratio as "X- 
ray faint." In these X-ray faint ellipticals, much of the 
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hot interstellar gas may have been lost in galactic winds 
(Loewenstein & Mathews 1987; David, Forman, & Jones 
1991), or by ram pressure stripping by ambient intraclus- 
ter or intragroup gas (White & Sarazin 1991). The source 
of the bulk of the X-ray emission in X-ray faint galaxies 
is uncertain; it might be due to low-mass X-ray binaries 
(LMXBs) like those seen in the bulge of our Galaxy (e.g., 
White, Nagase, & Parmar 1995), or to an active galactic 
nucleus (AGN) (Allen, di Matteo, & Fabian 2000), or to 
interstellar gas (Pellegrini & Fabbiano 1994), or to fainter 
stellar sources such as active M stars or RS CVn binaries 
(e.g., Pellegrini 1994). 

In general, X-ray faint galaxies exhibit significantly dif- 
ferent X-ray spectral properties than their X-ray bright 
counterparts. The X-ray bright galaxies are dominated by 
thermal emission at kT ~ 0.8 kcV due to their hot in- 
terstellar medium (ISM). On the other hand, the X-ray 
faint galaxies exhibit two distinct spectral components: 
first, they have a hard ~5-10 keV component, most eas- 
ily seen in ASCA spectra (Matsumoto et al. 1997), which 
is roughly proportional to the optical luminosity of the 
galaxy. Actually, both X-ray faint and X-ray bright early- 
type galaxies appear to have this hard X-ray component, 
which is roughly proportional to the optical luminosity of 
the galaxy. This suggests that the hard component is due 
to low-mass X-ray binaries (LMXBs) like those seen in 
the bulge of our Galaxy (e.g., White, Nagase, & Parmar 
1995). In some cases, an AGN may contribute to the hard 
component (Allen et al. 2000), or it might be due to hot 
gas. However, the ASCA observations do not resolve this 
component into discrete sources, nor to they provide much 
detailed information on its spectrum. 

X-ray faint galaxies also have a very soft (~0.2 keV) 
component, whose origin is uncertain (Fabbiano, Kim, & 
Trinchieri 1994; Pellegrini 1994; Kim et al. 1996). This 
difference in X-ray spectral characteristics is seen in the 
individual X-ray spectra of a number of X-ray bright and 
faint galaxies, and in the X-ray "colors" or hardness ratios 
determined for larger samples (Irwin & Sarazin 1998b). 
Suggested stellar sources for the very soft emission in X- 
ray faint ellipticals include active M stars, RS CVn bi- 
naries, or supersoft sources, but none of these appears to 
work quantitatively (Pellegrini & Fabbiano 1994; Irwin & 
Sarazin 1998a). It is possible that the soft X-rays are due 
to warm (0.2 keV) ISM (Pellegrini & Fabbiano 1994). Re- 
cently, we proposed that the very soft emission in X-ray 
faint ellipticals might be due to the same LMXBs responsi- 
ble for the hard emission (Irwin & Sarazin 1998a, b). Little 
is known about the very soft X-ray properties of Galactic 
LMXBs since most lie in directions of high Galactic hydro- 
gen column densities, so their soft X-ray emission is heavily 
absorbed. However, the origin of the very soft component 
in X-ray faint ellipticals remains something of a mystery. 

With the superb spatial resolution of the Chandra X-ray 
Observatory, it should be possible to resolve the emission 
from nearby X-ray faint early-type galaxies into LMXBs, 
if this is indeed the source of their emission. The clear test 
of the origin of the X-ray emission from X-ray faint ellip- 
ticals is to see if it resolves into discrete LMXBs. If the 
hard and/or soft X-ray emission from X-ray faint galaxies 
is due to LMXBs, it should come from a relatively small 
number of individually bright sources. Recently, we used 



a deep ROSAT HRI observation to detect a number of 
discrete X-ray sources in the nearby, optically luminous, 
X-ray faint elliptical NGC 4697 (Irwin, Sarazin, & Breg- 
man 2000, hereafter ISB). However, the bulk of the X-ray 
emission was not resolved. In ISB, we also simulated a 
40 ksec Chandra observation of NGC 4697, and showed 
that it should be possible to detect ^100 LMXBs if they 
provide the bulk of the emission. Here, we present the 
results of exactly this observation. At a distance of 15.9 
Mpc (Faber et. al. 1989; assuming a Hubble constant of 
50 km s" 1 Mpc -1 ), NGC 4697 is among the closest op- 
tically luminous, X-ray faint early-type galaxies. Given 
its proximity, NGC 4697 is an ideal target for detecting 
the LMXB population. It is sufficiently X-ray faint that 
diffuse ISM emission should not bury the emission of the 
fainter LMXBs. It should be possible to detect the LMXB 
population down to a level a luminosity of ~ 5 x 10 37 ergs 
s . The purposes of the Chandra observation are to re- 
solve and study the LMXB population of NGC 4697, to 
determine the source of both the hard and soft spectral 
components, and to detect or place strong limits on any 
residual diffuse (possibly gaseous) emission. 

Some initial results of this observations were presented 
in Sarazin, Irwin, & Brcgman (2000; hereafter Paper I). In 
the present paper, we give the detailed properties of the 
sources. In § 2, the observation and data analysis are de- 
scribed. The overall X-ray image is discussed in § 3. The 
properties of the resolved sources are given in § 4, while 
the remaining unresolved emission is analyzed in § 5. In 
§ 6, the X-ray spectral properties of the galaxy are derived. 
The interpretation of the results is discussed in § 7, and 
our conclusions are summarized in § 8. 

2. OBSERVATION AND DATA REDUCTION 

NGC 4697 was observed on 2000 January 15-16 with 
the ACIS-23678 chips operated at a temperature of -110 
C and with a frame time of 3.2 s. The pointing was deter- 
mined so that the entire galaxy was located on the S3 chip 
and so that the center of the galaxy was not on a node 
boundary of the chip. Although a number of serendipi- 
tous sources are seen on the other chips, the analysis of 
NGC 4697 in this paper will be based on data from the S3 
chip alone. The total exposure for the S3 chip was 39,434 
s. The data were telemetered in Faint mode, and only 
events with ASCA grades of 0,2,3,4, and 6 were included. 
We excluded bad pixels, bad columns, and the columns 
next to bad columns and to the chip node boundaries. We 
checked for periods of incorrect aspect solution, and none 
were found. The ACIS S3 chip is known to experience 
occasional periods of high background (Markevitch 2000a; 
Markevitch et al. 2000). We searched for such background 
flares by determining the lightcurve of the total count rate 
in the S3 chip binned in 40 s intervals; the total count rate 
is mainly due to background. No background flares were 
found, and the background rate was constant at about 0.9 
s _1 , which was the quiescent rate for the S3 chip at the 
time of the observation. Only events with photon energies 
in the range 0.3 to 10 keV were included in our analysis. 

This observation was processed at a time when the stan- 
dard pipeline processing introduced a boresight error of 
about 8" in the absolute positions of X-ray sources. We 
corrected for this using optical identifications and posi- 
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tions of X-ray sources. As noted below (§ 4), a number 
of the X-ray sources have faint optical identifications. Of 
these, three had accurate positions given in the USNO- 
A2.0 optical catalog (Monet et al. 1998). The optical and 
X-ray positions were offset by essentially the same amount 
and in the same direction (to within 0.3"). We applied the 
average of these offsets to the X-ray positions. We believe 
that the quoted absolute positions are accurate to about 
0.5" near the center of the S3 image, with larger errors 
further out. 

Determining the background for these observations 
proved to be difficult. We initially tried to use background 
from a series of nearly "blank sky" observations compiled 
by Markevitch (2000a,b). For the S3 chip, most of the 
exposure in these background files occurs at high Galactic 
latitude, mainly from the deep survey field of Mushotzky 
et al. (2000), These observations were taken when the fo- 
cal plane temperature was -110 C (as for our observations), 
and were screened to remove background flares; thus, they 
should be appropriate for our data. The blank-sky back- 
ground file gave a surface brightness in the hard (2-10 
keV) band which agreed well with the hard surface bright- 
ness in our data in outer parts of the S3 chip away from 
NGC 4697. However, our data had significantly more soft 
emission in the outer parts. In the outermost parts of the 
S3 chip, this extra soft emission had a fairly uniform sur- 
face brightness. While part of this emission is due to an ex- 
tended soft component of NGC 4697 (§ 5.3), the constant 
soft surface brightness component appeared to be an ad- 
ditional background component. We also determined the 
background for the region of the 13 chip which is furthest 
from NGC 4697 and along its minor axis. Our data also 
show a soft X-ray excess in this region relative to blank sky 
fields. We also examined the long ROSAT PSPC pointed 
observation of the galaxy (ISB), which showed excess soft 
emission in the outer parts of the field far from NGC 4697 
or other sources. 

On the sky, NGC 4697 is located near the edge of the 
North Polar Spur (NPS), a soft X-ray Galactic feature 
(e.g., Snowden et al. 1997). The ROSAT All Sky Sur- 
vey (RASS) images of this region (e.g., Snowden et al. 
1997) show a strong soft X-ray excess at the position of 
NGC 4697, which is consistent with the soft excess seen in 
the outer parts of the pointed PSPC observation. The ex- 
cess is most prominent in the R4R5 ROSAT band, which 
corresponds roughly to 0.4-1 keV. The excess in the Chan- 
dra data is mainly in our soft band (0.3-1 keV). The level 
of the soft excess due to the NPS is such that it would 
significantly affect the background in our observation. 

We used the surface brightness profile of NGC 4697 in 
various energy bands in both the ROSAT PSPC observa- 
tion and in our Chandra data to determine the portion 
of the extended soft emission which was due to the NPS 
and which was due to NGC 4697. Simple thermal mod- 
els with temperatures of ~ 2 x 10 6 K and a variety of 
abundances were used to convert the count rates for the 
soft excess between ROSAT and Chandra. In the outer 
parts of the S3 chip, we found that 55% of the soft excess 
was due to the NPS, while 45% was due to NGC 4697. 
For studying the properties of extended diffuse emission, 
we determine the background by combining the blank sky 
background of Markevitch (2000a, b) with background de- 



termined from the outer portion of the S3 chip in our ob- 
servation, and weighted by 45% and 55%, respectively. We 
added a systematic error to the background statistical er- 
ror based on the difference in these two determination of 
the background. 

In determining the properties of the resolved sources 
(§ 4), background for each source was determined from a 
region around the source. This background properly in- 
cludes all sources of diffuse emission, including the diffuse 
emission from NGC 4697. This also avoids the problem of 
the effect of the NPS as discussed above. 

All of the X-ray spectra were extracted using the PI 
values for the events in order to correct for gain varia- 
tions over the S3 chip. The PI values were recomputed 
using the gain file acisD1999-09-16gainN0004.fits, appro- 
priate for operating temperature of -110 C and frame- 
time of 3.2 s. The spectral responses were based on the 
FITS Embedded Function (FEF) files FP-110J31999-09- 
16fcf_piN0002. All of the spectra discussed here are ex- 
tracted from extended regions, which cover many of the 
32x32 pixel regions covered by individual FEFs. We used 
the RUNEXTRSPEC software package kindly provided by 
Alexey Vikhlinin (Vikhlinin, Markevitch, & Murray 2001), 
to extract spectra and to determine the response matrices. 
This program weights the response files by the number of 
counts in each region. 
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Fig. 1. — The Chandra S3 image of an approximately l'xl' re- 
gion near the center of NGC 4697. The greyscale varies with the 
square root of the X-ray surface brightness, which ranges from 1 to 
25 counts per pixel. (The ACIS pixels are 0'.'492 square.) The posi- 
tions of detected sources in the image are indicated by their source 
numbers in Table 1; the source numbers are ordered by increasing 
distance from the center of the galaxy. Src. 1 is located at the optical 
center of the galaxy. 

3. X-RAY IMAGE 

The raw Chandra S3 chip X-ray image was shown in 
Figure 1 of Paper I. This image showed the basic result 
of the Chandra observation: much of the emission from 
the galaxy is resolved into individual point sources of X- 
rays. The raw X-ray image of a 1' x 1' central region of the 
NGC 4697 is shown in Figure 1. The greyscale is propor- 
tional to the square root of the X-ray surface brightness, 
and the values vary between 1 and 25 cts pixel -1 . The 
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numbers identify the detected sources in this regions as 
listed in Table 1 below. The optical center of the galaxy 
coincides with the position of Src. 1. The center of the im- 
age has been offset from the optical because there are more 
sources to the ENE than WSW of the galaxy center. The 
density of resolved sources near the center of the galaxy is 
high, and they may become confused in this region. 




Fig. 2. — An adaptivcly smoothed version of the entire Chan- 
dra S3 X-ray image of the region around NGC 4697 in the 0.3—10 
keV band. The image was adaptively smoothed to a signal-to-noisc 
ratio of 3 per smoothing beam, and corrected for exposure and back- 
ground. Regions with an exposure of less than 20 ksec were removed 
to avoid artifacts at the chip edges. The grey scale is logarithmic, 
and ranges from 2 X 10~ 7 to 10~ 3 cts pixel -1 s _1 . (The ACIS pixels 
are 0'.'492 square.) 

In order to image the fainter, more diffuse emission, we 
adaptively smoothed the entire Chandra S3 X-ray image to 
a minimum signal-to-noise ratio of 3 per smoothing beam. 
The image was corrected for exposure and background. 
Regions with an exposure of less than 20 ksec were re- 
moved to avoid artifacts at the chip edges. This adaptively 
smoothed image is shown in Figure 2. The grey scale is 
logarithmic and covers surface brightnesses which range 
over a factor of >4000. The central 4'x4' region of the 
adaptively smoothed image was presented as Figure 2 in 
Paper I, showing only the brighter regions. 

For comparison, Figure 3 shows the Digital Sky Survey 
optical image of the same region as shown in Figure 2. The 
circles indicate the positions of individual X-ray sources 
which are identified in § 4 and Table 1 below. The dis- 
tribution of X-ray point sources in Figure 1 of Paper I 
appears to be elongated at the same position angle as the 
optical image of this E6 galaxy (67°; Jedrzejewski, Davies, 
& Illingworth 1987; Faber et al. 1989; Peletier et al. 1990); 
this is discussed in detail in § 4.6. (The small background 
galaxy LCRS B124537.4-053024, which is located at the 
right side of Figure 3, is beyond the edge of the S3 chip.) 

Figures 1, 2, and Figure 1 of Paper I show that the 
X-ray emission from NGC 4697 is mainly due to resolved 
sources, which are discussed in more detail in § 4 below. 
However, there is also some spatially extended, unresolved 
emission. Near the center of the image (Figure 2, Paper I), 
the unresolved emission follows the elliptical distribution 
of optical light in the galaxy. However, at larger radii the 



unresolved emission becomes somewhat irregular, and has 
something of an "L-shaped" distribution, with extension 
to the cast and south. The unresolved emission also is 
more extended than the optical light of the galaxy at the 
same relative surface brightness levels, particularly along 
the optical minor axis of the galaxy. The spatial distribu- 
tion of the unresolved emission is discussed in more detail 
in § 5.3, 
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Fig. 3. — The Digital Sky Survey optical image of NGC 4697, 
showing the same field as in Figure 2 or Figure 1 in Paper I. The 
circles show the positions of the discrete X-ray sources listed in Ta- 
ble 1. 

4. RESOLVED SOURCES 

4.1. Detections 

The discrete X-ray source population on the ACIS S3 
image (Figure 1 of Paper I) were determined, using us- 
ing both cell detection and wavelet detection algorithms. 
We used the CIAO CELLDETECT and wavdetect pro- 
grams. The two techniques gave consistent source lists 
with nearly identical fluxes. For sources located in the 
dense central region of NGC 4697 where confusion and 
the effect of an enhanced background due to unresolved 
sources or diffuse emission were important, the wavelet 
technique worked slightly better. The wavelet source de- 
tection threshold was set at 10~ 6 , which implies that <1 
false source (due to a statistical fluctuation in the back- 
ground) would be detected in the entire S3 image. This 
corresponded to requiring that the sources be detected at 
the >3-<7 level. The CELLDETECT threshold was initially 
taken to be >2-cr, but the resulting source list was culled 
to remove sources whose signal-to-ratio was less than 3- 
a. While variations in the background, exposure, and in- 
strumental Point-Spread-Function (PSF) across the image 
cause the minimum detectable flux to vary somewhat with 
position, over most of the galaxy the minimum detectable 
flux was about 2.6 x 10~ 4 cts s" 1 (L x = 5.0 x 10 37 ergs s" 1 
at the NGC 4697 distance) in the 0.3-10 keV band. Fluxes 
were corrected for exposure and the instrument PSF. The 
detection limit is higher at the edges of the field, and at the 
southern and western ends of the image, where the instru- 
mental PSF is larger than near the aim point. All of the 
source detections were verified by examining the image. 



CHANDRA X-RAY OBSERVATIONS OF NGC 4697 



5 



Table 1 
Discrete X-ray Sources 
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a The position of Src. 1 agrees with the optical center of NGC 4697 to within the errors. 
^Positions and fluxes of sources near the center of NGC 4697 are uncertain due to crowding. 

c This source appears to be extended, although this may be due to confusion with other sources or the effect of diffuse emission. 

d This source is at the edge of the S3 detector, and its flux is uncertain due to a large exposure correction. 

e Possiblc faint optical counterpart. 

■^Globular cluster is possible optical counterpart. 

g Source may be variable. 

^Supersoft source. 

lr Thc optical ID of this source is an AGN at rcdshift z — 0.696. 
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Table 1 lists the 90 discrete sources detected by this 
technique, sorted in order of increasing distance d from the 
center of NGC 4697. Columns 1-9 give the source number, 
the IAU name, the source position (R.A. and Dec, J2000), 
the projected distance d from the center of NGC 4697, the 
projected isophotal semi-major axis a (see § 4.6 below), 
the count rate and the \-a error, the signal-to-noise ratio 
SNR for the count rate, and the unabsorbed 0.3-10 keV 
X-ray luminosity Lx (assuming the source is located at 
the distance of NGC 4697, see § 4.3 below). The conver- 
sion to luminosities assumes the best-fit source spectrum 
(thermal bremsstrahlung with kT — 8.1 keV and Galactic 
absorption with N H = 2.12 x 10 20 cirT 2 ; § 6 and Table 4, 
row 3). The statistical errors in the positions of most of 
the sources are quite small (<~ 0.2"), and the overall abso- 
lute errors are probably ^0.5" near the center of the field, 
with larger errors near the outside of the field. The val- 
ues of the distance from the galaxy center d are actually 
computed from the position of Src. 1; its position agrees 
with the position of the optical center of the NGC 4697 
to within the errors, but relative X-ray positions are more 
accurate. Thus, d might differ from the distance from the 
optical center of the galaxy by <~1". 

As noted above, our detection limit for sources should 
result in <1 false source (due to a statistical fluctuation) in 
the entire S3 field of view. However, many of the detected 
sources may be unrelated foreground or (more likely) back- 
ground objects. Based on the source counts in Brandt 
et al. (2000) and Mushotzky et al. (2000), we would ex- 
pect about 10-15 serendipitous sources in our observation. 
These should be spread out fairly uniformly over the S3 
image (Figure 1 of Paper I), except for the effects of the 
reduced exposure and increased PSF at the outside of the 
field. Thus, the unrelated sources should mainly be found 
at larger distances from the optical center of NGC 4697 
(the bottom part of Table 1), while the sources associated 
with NGC 4697 should be concentrated to the center of 
the galaxy. Within 2' of the center of NGC 4697 (roughly 
the region covered by Figure 2 in Paper I), ~2 of the ~60 
the detected sources would be expected to be unrelated to 
NGC 4697. 

4.2. Identifications 

The position of Src. 1 agrees with the optical position 
of the center of NGC 4697 (R.A.=12 h 48 m 35!71, Dec .=- 
5°48'02'.'9; Wegner et al. [1996]) to within the combined 
X-ray and optical errors. Its luminosity is Lx — 8 x 10 38 
ergs s -1 . This source appears to be broader than the in- 
strumental PSF (Fig. 1). On the other hand, the density 
of sources is quite high near the center of NGC 4697, and 
the apparent extension may be due to confusion with other 
sources or diffuse emission. It is possible that the central 
source is due in part to an AGN and/or to one or more 
LMXBs, as discussed below in § 7.1. 

Near the center of the galaxy, the density of sources 
is quite high, and it is likely that some of the positions, 
fluxes, and sizes are affected by source confusion. Src. 3 is 
very close to Src. 1, and also appears somewhat extended. 
It may be a composite or may be affected by the central 
Src. 1. Src. 2 is very close to both Src. 1 and Src. 4. Srcs. 7 
and 10 are quite close together. 

ISB detected 12 X-ray sources around NGC 4697 using 



the ROSATEKL, all of which lie within the field of view of 
the S3 image. All but one of these sources lie within a few 
arcsec of strong sources which we detect with Chandra. 
The identifications of these sources (the sources numbers 
in Table 1 from ISB) are listed in column 10 of Table 1. 
ISB Src. 7 corresponds to the central X-ray source (Src. 1 
in the present paper). ISB Src. 8 was an extended source 
located about 7" north of the center of the galaxy; we be- 
lieve this source is actually due to the superposition of our 
sources 6 and 8, with possible contributions from sources 
5 and/or 9 as well. In general, the ISB sources have Chan- 
dra S3 count rates > 2 x 10~ 3 cts s -1 , and are thus almost 
an order of magnitude brighter than the faintest sources 
we detect. ISB Src. 3 was located about 54" west of the 
center of NGC 4697. We don't detect a source or any en- 
hancement in the background surface brightness at this 
location. Given that ISB Src. 3 was the least significant 
source detected by ISB (only slightly better than 2-cr), it 
seems likely that this source was a statistical fluctuation in 
the ROSATBBI image. Alternatively, this source may be 
variable, and might have declined by more than an order 
of magnitude since the ROSAT observation (§ 4.5). 

The positions of the X-ray sources were examined on the 
Digital Sky Survey (DSS) image of this region (Fig. 3). 
Eleven of the sources had possible faint optical counter- 
parts on this image; four others appear on deeper optical 
images discussed below. These are all marked with note 
"e" in Table 1. The possible faint optical counterparts of 
Srcs. 71 and 72 were noted previously by ISB. None of the 
sources corresponded to objects list on NED or SIMBAD. 

However, several of the X-ray sources coincide with can- 
didate globular clusters associated with NGC 4697. Src. 64 
has a faint optical counterpart on the DSS image which 
agrees with the position of globular cluster candidate 33 
from the catalog of Hanes (1977). (However, this source 
has an unusual X-ray spectrum, which may indicate that 
is is actually a background AGN; see § 4.4 below.) Most 
of the Hanes (1977) globulars are fairly bright and located 
at relatively large distances from NGC 4697. J. Kave- 
laars (2000, private communication) gives a list of gener- 
ally fainter candidate globular clusters located between 1.'5 
and 2.'5 from the center of NGC 4697. By comparing this 
list with the X-ray source positions we find faint candidate 
globular clusters located within 1" of the positions of X-ray 
Srcs. 53, 54, 55, 61, 64, 67, and 69. The globular cluster 
associated with Src. 64 is the same one listed by Hanes 
(1977). Given the density of candidate globulars and the 
very good position agreement, one would expect ^0.3 as- 
sociations to occur at random. Thus, it is likely that all of 
the 7 associations of X-ray sources with optical globular 
cluster candidates in Table 1 are real. However, at the dis- 
tance to NGC 4697, globular clusters are not resolved in 
ground-based optical images, and the candidate globulars 
were identified by luminosities and (possibly) colors. As a 
result, as many as half of them might be unrelated faint 
optical objects, rather than globular clusters. 

An example of such a miss-identified globular cluster 
is the the optical object associated with Src. 72. Hanes 
(1977) lists this object as a candidate globular cluster 
(# 24), although it would be rather bright for a globu- 
lar at the distance of NGC 4697. Moreover, the USNO- 
A2.0 optical catalog (Monet et al. 1998) indicates that this 
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Table 2 
Hardness Ratios for Sources 



Src. 




H21 




H31 


Src. 




H21 




H31 


Src. 




H21 




H31 


1 


-0.22 


-0.31.-0.13) 


-0.34 


-0.42,-0.25) 


32 


+ 0.02 


-0.17, + 0.21) 


-0.33 


-0.52,-0.15) 


62 


+ 0.00 


-0.30. +0.30) 


-0.21 


-0.51, + 0.10) 


2 


+ 0.01 


-0.32, +0.34) 


-0.68 


-1.00,-0.36) 
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-0.30 


-0.70, +0.10) 
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-1.00,-0.22) 


63 


-0.38 


-0.65,-0.12) 
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-0.46, + 0.10) 
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-0.20, +0.10) 


-0.92 


-1.00,-0.78) 


4 


-0.03 


-0.33, + 0.27) 


-0.42 


-0.71,-0.14) 


35 


-0.18 


-0.39, +0.04) 


-0.53 


-0.75,-0.31) 


65 


+0.24 


+0.14, + 0.35) 


-0.30 


-0.44,-0.15) 


5 


-0.21 


-0.55, + 0.12) 


-0.09 


-0.44, +0.26) 


36 


+0.03 


-0.06, + 0.12) 


-0.26 


-0.35,-0.18) 


66 


+ 0.07 


-0.17. +0.31) 


-0.83 


-1.00.-0.55) 


6 


4-0.00 


-0.16, +0.16) 


-0.32 


-0.47,-0.18) 


37 


-0.16 


-0.37, + 0.06) 


-0.67 


-0.90,-0.45) 


67 


+ 0.10 


-0.87, + 1.00) 


+ 0.39 


-0.23, + 1.00) 


7 & 10 a 


-0.11 


-0.33, + 0.12) 


-0.67 


-0.90,-0.44) 


38 


-0.21 


-0.55, +0.12) 


-0.35 


-0.70, +0.00) 


68 


+ 0.30 


-0.12, + 0.71) 


-1.00 


-1.00,-0.07) 


8 


+ 0.17 


+0.01, + 0.32) 


-0.01 


-0.21, +0.19) 


39 


+ 0.19 


-0.07, + 0.45) 


-0.25 


-0.60, + 0.10) 


69 


+ 0.07 


-0.07, + 0.21) 


-0.22 


-0.37,-0.06) 


9 


-0.62 


-0.74,-0.49) 


-0.87 


-1.00,-0.75) 


40 


+0.02 


-0.07, + 0.12) 


-0.23 


-0.32,-0.14) 


70 


-0.10 


-0.56, +0.36) 


-0.04 


-0.53, + 0.45) 


11 


-0.07 


-0.40, + 0.27) 


-0.34 


-0.68, +0.00) 


41 


-0.14 


-0.32, +0.04) 


-0.47 


-0.65,-0.29) 


71 


-0.23 


-0.36,-0.10) 


-0.67 


-0.81.-0.54) 


12 


-0.17 


-0.42, + 0.07) 


-0.45 


-0.70,-0.20) 


42 


+0.00 


-0.18, + 0.18) 


-0.21 


-0.38,-0.04) 


72 


+0.07 


+0.01, + 0.14) 


-0.31 


-0.37,-0.24) 


13 


-0.27 


-0.57, + 0.04) 


-0.98 


-1.00,-0.59) 


43 


-0.23 


-0.47, +0.01) 


-0.25 


-0.50, +0.00) 


73 


-0.38 


-0.65,-0.12) 


-0.69 


-0.99,-0.38) 


14 


-0.08 


-0.23, + 0.06) 


-0.49 


-0.62,-0.35) 


44 


-0.35 


-0.55,-0.15) 


-0.22 


-0.43,-0.01) 


74 


-0.27 


-0.46,-0.07) 


-1.00 


-1.00,-0.80) 


15 


— 0.50 


— 0.86,-0.14) 


— 0.29 


— 0.68, +0.09) 


45 


— 0.12 


— 0.33, +0.10) 


— 0.08 


— 0.31, + 0. 14) 


75 


— 0.37 


— 0.59, — 0. 16) 


— 0.66 


— 0.90,-0.43) 


16 


-1.00 


-1.00,-0.94) 


-1.00 


-1.00,-0.93) 


46 


+0.00 


-0.87, + 0.87) 


+ 0.24 


-0.36, +0.85) 


76 


-0.04 


-0.23, +0.16) 


-0.19 


-0.38, + 0.00) 


17 


+ 0.00 


-0.32, +0.32) 


-0.22 


-0.54, +0.10) 


47 


-0.06 


-0.26, +0.15) 


-0.14 


-0.35, +0.06) 


77 


-0.42 


-0.51,-0.34) 


-0.63 


-0.72.-0.55) 


18 


-0.46 


-0.78,-0.14) 


-0.25 


-0.60, +0.09) 


48 


+0.10 


-0.37, + 0.57) 


-0.05 


-0.62, +0.52) 


78 


-0.06 


-0.20, +0.07) 


-0.55 


-0.68,-0.42) 


19 


-0.95 


-1.00,-0.89) 


-1.00 


-1.00,-0.96) 


49 


-0.27 


-0.57, + 0.04) 


-0.64 


-0.99,-0.30) 


79 


+ 0.00 


-0.53, +0.52) 


+ 0.11 


-0.36, +0.57) 


20 


+0.36 


+0.18, + 0.55) 


-0.08 


-0.40, +0.25) 


50 


-0.01 


-0.38, +0.36) 


-0.92 


-1.00,-0.47) 


80 


+ 0.37 


-0.10, +0.85) 


-0.09 


-0.99, + 0.82) 


21 


-0.62 


-0.90,-0.33) 


-0.44 


-0.74,-0.14) 


51 


+0.15 


-0.22, + 0.52) 


+ 0.13 


-0.27, +0.52) 


81 


+0.50 


+ 0.21, + 0.79) 


+ 0.24 


-0.21, +0.70) 


22 


-0.05 


-0.21, +0.11) 


-0.05 


-0.21, +0.12) 


52 


-1.00 


-1.00,-0.94) 


-1.00 


-1.00,-0.92) 


82 


+0.82 


+0.79, + 0.85) 


+0.83 


+0.80, + 0.86) 


23 


-0.53 


-0.80,-0.27) 


-0.24 


-0.52, +0.04) 


53 


-0.10 


-0.21, +0.02) 


-0.48 


-0.59,-0.37) 


83 


+ 1.00 


-1.00, + 1.00) 


+ 1.00 


+0.68, + 1.00) 


24 


-0.21 


-0.35,-0.07) 


-0.48 


-0.61,-0.34) 


54 


-0.41 


-0.77,-0.05) 


-0.42 


-0.81,-0.02) 


84 


+ 0.15 


-0.14, +0.45) 


-0.05 


-0.44, + 0.34) 


25 


-0.44 


-0.70,-0.17) 


-0.31 


-0.59,-0.03) 


55 


+0.48 


-0.09, + 1.00) 


+ 0.45 


-0.16, + 1.00) 


85 


+ 0.05 


-0.24, +0.33) 


-0.98 


-1.00.-0.63) 


26 


-0.04 


-0.31, + 0.22) 


-0.28 


-0.54,-0.01) 


56 


+0.00 


-0.55, + 0.54) 


-0.05 


-0.61, +0.52) 


86 


-0.35 


-0.74, +0.04) 


-0.69 


-1.00,-0.20) 


27 


+ 0.17 


+0.07, + 0.28) 


-0.08 


-0.21, +0.06) 


57 


+0.03 


-0.19, + 0.25) 


-0.26 


-0.49,-0.03) 


87 


-0.74 


-1.00, +0.15) 


+ 0.29 


-0.24, + 0.82) 


28 


-0.26 


-0.45,-0.07) 


-0.47 


-0.66,-0.27) 


58 


+ 0.35 


+ 0.04, +0.67) 


-0.04 


-0.61, +0.52) 


88 


+0.53 


+0.19, + 0.87) 


+ 0.20 


-0.41, + 0.80) 


29 


-0.28 


-0.55,-0.01) 


-0.69 


-0.99,-0.40) 


59 


+0.00 


-0.42, + 0.42) 


-0.71 


-1.00,-0.23) 


89 


+ 0.04 


-0.20, +0.29) 


-0.89 


-1.00.-0.55) 


30 


+ 0.00 


-0.29, +0.28) 


-0.30 


-0.58,-0.03) 


60 


-0.12 


-0.66, +0.41) 


-0.05 


-0.62, +0.52) 


90 


+0.89 


+0.65, + 1.00) 


+0.86 


+0.54, + 1.00) 


31 


+ 0.13 


-0.21. +0.47) 


-0.12 


-0.54, +0.30) 


61 


-0.44 


-0.58,-0.31) 


-0.56 


-0.70,-0.43) 













a Srcs. 7 & 10 are very close together, and we report their combined hardness ratios. 



source has an very blue color (B — R « —0.7), which is 
bluer than even a Raylcigh-Jeans spectrum, and would be 
unheard of for a globular cluster. An optical spectrum 
kindly provided by C. Mullis (2000, private communica- 
tion) shows that this object is a background AGN at a 
rcdshift of z = 0.696, rather than a globular cluster. The 
very blue color is due in part to a strong Mg II A2798 line 
in the spectrum, rcdshiftcd into the B band. 

4.3. Luminosities and Luminosity Function 

The count rates for the sources were converted into un- 
absorbed luminosities (0.3-10 keV) assuming that all of the 
sources were at the distance of NGC 4697, which we take 
to be 15.9 Mpc (Faber et al. [1989] and H Q = 50 km s" 1 
Mpc -1 ). We adopted the best-fit Chandra X-ray spectrum 
of the resolved sources within the inner 1 effective radius 
(Table 4, row 3 below). The resulting factor for convert- 
ing the count rate (0.3-10 keV) into the unabsorbed lumi- 
nosity L x (0.3-10 keV) was 1.92 x 10 41 ergs cnt" 1 . The 
resulting X-ray luminosities are given in column 9 of Ta- 
ble 1, and range from about 5 x 10 37 to 2.5 x 10 39 ergs s _1 . 
The cumulative luminosity function of all of the sources 
was shown as a histogram in Figure 3 of Paper I. There 
we found that the luminosity function of the NGC 4697 
sources could not be fit by a single power law, as the lu- 
minosity function has a "knee" at Lx ~ 3 x 10 38 ergs s _1 . 
A broken power-law, 

dN { (L 38 /L b )- a < L 38 < L b 

Tj~ =N ° { > (!) 

dL38 { (L 38 /L b )-^ L 38 >L b 

gave a good fit, where L 38 is the X-ray luminosity (0.3-10 
keV) in units of 10 38 ergs s _1 . The best fit, determined 
by the maximum-likelihood method, gave N Q — 8.OI5 2, 
a t = 1.29lo; 3 g, a h = 2.76tJ; 8 g, and a break luminosity 
of L b = 3.2to° 8 x 10 38 ergs s^ 1 We also determined the 



luminosity functions for sources within the optical half- 
light elliptical isophote and for sources are larger distances. 
The errors were increased significantly (particularly for the 
outer region), but the luminosity functions agreed with the 
one shown in Figure 3 of Paper I. 

4.4. Hardness Ratios of Sources 

We studied the crude spectral properties of the resolved 
sources by using hardness ratios. Hardness ratios or X- 
ray colors have the advantage that they can be applied to 
weaker sources. We defined two hardness ratios as H21 
= {M-S)/{M + S) and H31 = (H — S)/(H + S), where 
S, M, and H are the total counts in our soft (0.3-1 keV), 
medium (1-2 keV), and hard (2-10 keV) bands, respec- 
tively. The hardness ratios of the sources and 90% con- 
fidence regions are listed in Table 2 (see also Figure 4 in 
Paper I). For comparison, the hardness ratios (H21,H31) 
are (-0.38, -0.57) for all of the emission, (-0.69, -0.82) 
for the unresolved emission, and (—0.14, —0.37) for the 
sum of the sources, all within one effective radius (§ 5). 

There are three moderate luminosity sources with hard- 
ness ratios of (—1, —1), which means that they have no de- 
tectable emission beyond 1 keV. For example, a blackbody 
spectrum with a temperature of 0.1 keV and the Galactic 
absorption towards NGC 4697 (N H = 2.12 x 10 20 cur 2 ) 
would give hardness ratios of (H21,H31) w (-0.97, -1.00). 
These three sources, Src. 16, Src. 19, and Src. 52, are al- 
most certainly supersoft sources (e.g., Kahabka & van den 
Heuvel 1997). All three are located at small enough dis- 
tances from the center of NGC 4697 that they are likely to 
be associated with the galaxy. Srcs. 19 & 52 also may be 
variable (§ 4.5). We believe this observation represents the 
first detection of supersoft sources in a luminous elliptical 
galaxy. 

There are three sources with hardness ratios of ~(1, 1), 
one of which is the brightest source in the field (Src. 82), 
These sources are Src. 82, Src. 83, and Src. 90. As is clear 
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from their source numbers, all of these sources are located 
far (>3.'5) from the center of NGC 4697. These are prob- 
ably unrelated, strongly absorbed AGNs, similar to the 
sources which produce the hard component of the X-ray 
background, and which appear strongly at the faint fluxes 
in the deep Chandra observations of blank fields (Brandt et 
al. 2000; Mushotzky et al. 2000; Giacconi ct al. 2001). As a 
comparison, a power-law spectrum with a photon spectral 
index T = 1.5 and with an absorbing column of Nh = 10 22 
cm" 2 gives hardness ratios of (H21,H31) « (+0.88, +0.92). 

There are eight sources (Srcs. 13, 50, 64, 66, 68, 74, 85, 
& 89) with hardness ratios near (0, —1) which have essen- 
tially no hard emission. For example, a power-law spec- 
trum with a photon spectral index T = 3.5 and with an 
absorbing column of Njj = 3 x 10 21 cm" 2 gives hardness 
ratios of (H21,H31) « (-0.18, -0.81). Six of these sources 
are at large radii (>2.'2), which suggests that most of this 
population is also unrelated to NGC 4697. The hardness 
ratios for these sources are consistent with the soft-band- 
only sources seen in the Giacconi et al. (2001) deep field 
at similar flux levels. However, note that Src. 64 was also 
identified with a candidate globular cluster in NGC 4697; 
perhaps this source is a globular cluster LMXB with an 
unusual X-ray spectrum, or the optical counterpart may 
be a background AGN rather than a globular. (Note that 
the 11 possible background sources selected by colors in 
the last two paragraphs plus the AGN [Src. 72] would ac- 
count for most of the ~10-15 unrelated sources expected 
based on deep blank sky observations.) 

Most of the sources lie in a diagonal swath centered 
at about (—0.15,-0.40). Note that hardness ratios of 
(—0.14, —0.37) correspond to the cumulative source spec- 
trum with one effective radius (Table 4, row 3 below). 
This is a hard spectrum fit by thermal bremsstrahlung 
with kT — 8.1 keV and Galactic absorption with Nh = 
2.12 x 10 20 cm" 2 . These values are similar to but slightly 
harder than the integrated colors or the entire galaxy, but 
considerably harder than the values for the unresolved 
emission. There may be some tendency for the fainter 
sources to have higher H/M ratios than the brighter 
sources. 

4.5. Variability of Sources 

We searched for variability in the X-ray emission of the 
resolved sources over the duration of the Chandra obser- 
vation using the KS test. This test can detect a variation 
in the flux from the source over the ~11 hour duration of 
the observation, such as a secular increase in the flux or a 
sudden turn-on or off of the source. In terms of the phe- 
nomenology of LMXBs in our Galaxy, one might detect 
orbital variations or other secular variations. Although we 
could, in principal, detect a Type I X-ray burster, the lu- 
minosities (typically < the Eddington luminosity of a NS) 
and durations (10-10 3 s) imply that we would expect < 
1 count from a burst at the distance of NGC 4697. The 
KS test cannot detect short term variations in the source 
which do not individually contribute significantly to the 
total number of counts. For example, we could not detect 
periodic pulses due to rotational modulation of an accret- 
ing neutron star. In principal, other tests could be done 
to search for pulsations from some of the brighter sources; 
in practice this is difficult because of the 3.2 s periodic- 



ity on the events imposed by the ccd readout. We found 
that 11 sources were inconsistent with a constant flux at 
>90% confidence; these sources are noted in Table 1. In 
Figure 4, we show the histogram of the total counts from 
the five sources with more than 20 net counts and with 
variability detected at >95% confidence. 
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Fig. 4. — The solid histograms gives the accumulated fraction of 
events for sources as a function of the accumulated exposure time. 
The dashed line is the predicted distribution under the hypothesis 
that the source plus background rate is constant. The five sources 
shown are those with more than 20 net counts (subtracting back- 
ground) and with variability detected at >95% confidence. 

We also tested for variability over longer time scales by 
comparing the Chandra count rates of the sources with 
detections or limits from a long (78,744 s) ROSAT HRI 
observation of the galaxy on 1997 June 19 - July 20 (ISB). 
We considered the sources as variable if they increased 
or decreased by >50% at >lcr significance. Sources with 
significantly higher fluxes at the time of the Chandra ob- 
servation were Srcs. 19, 27, 52, 61, 64, 72 (ISB Src. 11), 78, 
and 82. Src. 64 also appeared to vary during the Chandra 
observation. Both Srcs. 19 and 52 are supersoft sources. 
The third supersoft source (Src. 16) also showed some ev- 
idence for variability, but not at as significant a level. 

4.6. Spatial Distribution of Sources 

One might expect the stellar X-ray sources associated 
with NGC 4697 to have a spatial distribution which is 
very similar to that of the optical light. Unfortunately, the 
spatial distribution of detected X-ray sources may also be 
affected by variations in the sensitivity limit due to vari- 
ations in exposure, the instrumental PSF, and the diffuse 
background and emission. Because of these possible prob- 
lems and statistical limits due to the number of sources, 
we limit this discussion to a few simple comparisons of the 
X-ray source and optical distributions. We consider only 
sources within 3'.5 of the center of NGC 4697, which is 
slightly less than the distance to the nearest chip edge. 
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Figure 5 shows the distribution of the absolute values of 
the position angle (PA) of the sources. Here, PA is mea- 
sured from north to east, and is taken to be in the range 
— 180° < PA < 180°. For a set of sources whose projected 
density is constant on concentric aligned ellipses of con- 
stant ellipticity e, the distribution of angles <f> relative to 
the semimajor axis is 

JV(1 - e) d<t> 



n{4>) d(f> 



2ir (1 - e )2 



(2) 



where n(<p) d<f> is the number of sources with angles be- 
tween <f> and 4> + d<f>, and N is the total number of sources. 
The short-dash curve in Figure 5 is the predicted distri- 
bution based on the optical photometry of the galaxy, 
plus the expected number of background sources in this 
region. The optical photometry gives PA — 67° ± 4° and 
e = 0.42 ±0.06 (Jedrzejewski et al. 1987; Faber et al. 1989; 
Pelcticr et al. 1990); the errors approximately cover the 
range of values in the literature plus the radial variation 
at radii greater than 10". The long-dash curve is the best- 
fit elliptical distribution plus background, determined by 
a maximum-likelihood fit to the observed values. The best 
fit distribution gives PA = 79° ± 15° and e = 0.40 ± 0.09. 
The best fit curve is not a statistically significant improve- 
ment on the optical model. 



did this analysis using the elliptical semimajor axes of the 
sources (a values in Table 1), but the results were nearly 
identical. These comparisons show that the X-ray source 
distribution is elongated in about the same direction and 
by about the same amount as the optical light, and that 
the radial distribution of X-ray sources is proportional to 
that of the optical light, all to within the errors. 



1.2 1.4 1.6 
log r (arcsec) 



1.8 



Fig. 6. — The solid histogram gives the accumulated number of 
X-ray sources within 3.'5 of the center of NGC 4697 as a function 
of the radius. The dash curve is the predicted distribution based on 
the optical photometry of the galaxy, plus the expected number of 
background sources in this region. 
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Fig. 5. — The solid histogram gives the distribution of the posi- 
tion angles of the X-ray sources within 3f5 of the center of NGC 4697 
in 10° bins, as a function of the absolute value of the position angle 
PA. PA is measured from north to east, and is taken to be in the 
range -180° < PA < 180°. The short-dash curve is the predicted 
distribution based on the optical photometry of the galaxy, plus the 
expected number of background sources in this region. The long- 
dash curve is the best-fit elliptical distribution plus background, 
determined by a maximum-likelihood fit to the observed values. 

We also compared the radial distribution of sources with 
the photometry of the optical light. Figure 6 shows the ac- 
cumulated source number as function of the radius. The 
dashed curve is the predicted distribution if the source 
counts follow the optical light in the galaxy plus the ex- 
pected number of background sources. The optical pho- 
tometry was modeled as a de Vaucouleurs profile with ef- 
fective semimajor and semiminor axes of a or f = 95" and 
b e g = 55". The KS test indicates that the optical distri- 
bution is a good fit to the X-ray source distribution. We 
determined the best-fit de Vaucouleurs' profile plus back- 
ground fit to the source radial distribution, but is was 
not distinguishable from the optical distribution. We also 



5. DIFFUSE EMISSION 

5.1. Resolved vs. Diffuse Emission 

We determined the portion of the X-ray emission due 
to resolved sources and unresolved emission, both for the 
entire Chandra X-ray band (0.3-10 kcV) and for three nar- 
rower bands: hard H (2-10 keV), medium M (1-2 keV), and 
soft S (0.3-1 keV). We found that these bands gave reason- 
able count rates and spectral discrimination for a spectrum 
like that of the entire galaxy (§ 6), which has very soft and 
very hard components. The emission was determined for 
two spatial regions. First, it was done within the ellipti- 
cal optical isophote which contains one half of the optical 
light from the galaxy. We refer to these counts as coming 
from within "one effective radius." This isophote has an 
effective semimajor axis of 95", an effective semiminor axis 
of 55", an ellipticity of 0.42, and a position angle of 67°. 
We also determined the fluxes within an elliptical annular 
region ranging from one to two effective radii. The popula- 
tion of resolved sources was determined as discussed above 
(§ 4). Then, the total emission was determined, and the 
source flux subtracted to give the amount of unresolved 
emission. 

The results are listed in Table 3. For each spatial re- 
gion, the emission is divided into "Resolved" and "Unre- 
solved" components, with "AH" denoting the sum of these 
two. The truly "Diffuse" emission is the unresolved emis- 
sion, corrected for unresolved binary sources (§ 5.2), while 
"Discrete" emission is the sum of resolved and unresolved 
binary sources. The emission is divided into the soft (S), 
medium (M), and hard (H) bands, with "Total" denoting 
the sum of these (all emission from 0.3 to 10 keV). In ad- 
dition to the total number of counts, we give the hardness 
ratios H21 and H31 as defined in S 4.4. The fractional con- 
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Table 3 
Resolved vs. Diffuse Emission 









Hardness 


Fraction 


Luminosity 


Region Origin 


Band 


Counts 


(H21 or H31) 


(%) 


(10 40 ergs s" 1 ) 


< 1 a e ff All 


Total 


3715+90 




(100) 


1.57 




H 


584+69 


-0.57+°°4 

— U.U4 


(100) 






M 


975+36 


-0.38t°;°g 


(100) 






s 


2154±107 




(100) 




J. l^^OVJl v \_>v_x 


Total 


2271±57 




61+10 


1.11 




H 


474+31 


-0 37+2 ^ 

" J ' -0.03 


81+14 






M 


771 -J-Qc: 


-0 14+ 003 








c 
o 




/1 8-l-Q 




TT 

U nresolved 


rorai 


1 /l/l/l-t-1 07 






U.4D 




jj 


1 1 n+7fi 

± ±uzn ( u 




1 Q+1 4 






M 




-0 69+ ' 08 


21+5 






Q 


1 1 98+1 1 ^ 


^9+^ 




JJlIIUSC 


rorai 

s 


OCQ 1 1 QQ 

869+1 38 

WU —1 L tJU 




9Q-I-/1 
40+7 


U. lo 


Discrete 


Total 


2846+166 




77±4 


1.39 
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tribution of resolved, unresolved, and truly diffuse emis- 
sion to the counts within each band is given. Finally, the 
unabsorbed luminosity for each component and band are 
given, based on the best-fit spectra for each component 
(§ 6 and Table 4 below). 

Within one effective radius, 61% of the counts and 71% 
of the X-ray luminosity is resolved into individual X-ray 
sources for the total band. In the hard, medium, and soft 
bands, the resolved count fractions are 81%, 79%, and 
48%, respectively. Between 1 and 2 effective radii, the re- 
solved fraction drops but the spectrum of the unresolved 
emission gets harder, suggesting that this is partly the re- 
sult of lower sensitivity to point sources at radii where the 
point spread function (PSF) is larger. However, it is also 
clear that there is a very spatially extended soft compo- 
nent to the X-ray emission. 

5.2. Unresolved Binary Sources and Truly Diffuse 
Emission 



A portion of the unresolved emission must also come 
from LMXBs which fall below our threshold for reliable 
detection of resolved sources. Indeed, if one sets the de- 
tection threshold for sources lower, one finds many more 
2-<7 fluctuations than expected just from Poisson statis- 
tics. However, it would be difficult to use such a fluctua- 
tion analysis to extend the statistical detection of LMXBs 
down a factor of ~50 to Lx ~ 10 36 ergs s _1 , the approx- 
imate lower limit for LMXBs in globular clusters (Hertz 
& Grindlay 1983). If the observed luminosity function of 
the resolved sources in NGC 4697 (Fig. 4, Paper I) is ex- 
tended down to Lx = 10 36 ergs s , the contribution of 
LMXBs to the total band X-ray emission from NGC 4697 
increases by about 17%. This result agrees approximately 
with the result from the X-ray colors of the sources and 
unresolved emission given in the next paragraph. Still, 
this is a rather significant extrapolation, which could be 
inaccurate if the slope of the luminosity function changes 
below Lx = 5 x 10 37 ergs s _1 , as is seen in the bulge of 
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M31 (Shirey et al. 2001). 

We have estimated the portion of the unresolved emis- 
sion which is due to unresolved LMXBs based on the X- 
ray hardness ratios of the sources. The spectral analy- 
sis (§ 6) indicates that LMXBs produce the majority of 
the X-ray luminosity at photon energies above 1 kcV. 
Thus, we assume that all of the hard (H) and medium 
(M) counts in the unresolved emission come from unre- 
solved discrete sources. Moreover, we assume that these 
unresolved sources have the same spectral properties as 
the resolved sources. We use the hardness ratios of the 
resolved sources and the H+M counts of the unresolved 
emission to correct the soft (S) band unresolved emission 
for unresolved LMXBs. The counts and fractions for this 
truly "Diffuse" component are listed in Table 3. When 
corrected in this way, it appears that LMXBs (resolved 
and unresolved) provide 77% of the total band counts and 
89% of the luminosity in the inner effective radius. 

The X-ray spectrum of the unresolved emission consists 
of two components: a hard component with a spectrum 
which is consistent with that of the resolved sources, and 
a soft component (§ 6.2 below). If we attribute the hard 
component to unresolved sources, then we can determine 
the total luminosity of the discrete sources from the lumi- 
nosity of the hard component in the total spectrum (§ 6 
below). Based on the spectral fits, we would conclude 
that 74% of the total luminosity is due to discrete sources. 
The differences in the results using colors or spectral fits 
partly reflects uncertainties in the spectral fits and count 
rates, and differences in the spectral models of the differ- 
ent components. However, this difference may also be due 
to errors in the calibration of the ACIS S3 below 0.7 keV 
(§ 6 below). As a result of concerns about the low en- 
ergy spectral calibration, we use the hardness ratios and 
counts rather than spectra to determine the portion of 
the unresolved emission which is due to LMXBs. Using 
the spectra requires extrapolating the fit for the spectral 
range 0.7-10 keV down to 0.3 keV, where the spectral cal- 
ibration appears to be very uncertain. Using the hardness 
ratios involves only assuming that the unresolved LMXBs 
have the same spectrum as resolved sources. 

The remaining ~11% of the luminosity and ^23% of the 
counts would come from a more diffuse component with 
a soft (~0.3 keV) spectrum. In § 7.3, we show that this 
emission is almost certainly due to diffuse gas. 

5.3. Spatial Distribution of Diffuse Emission 

We determined the radial distribution of the diffuse 
emission. Because we are mainly interested in the (pre- 
sumably gaseous) soft component, we did this in the soft 
band (0.3-1.0 keV). Resolved sources were excluded. The 
unresolved emission was collected in circular annuli. (We 
also determined the surface brightness in elliptical annuli 
whose orientation and ellipticity matched that of the opti- 
cal isophotes of the galaxy, and the results were essentially 
identical to those for circular annuli.) The outer bound- 
ary of the largest annulus was the largest circle which fit 
entirely on the S3 chip; thus, regions of the chip corners, 
south edge, and west edge were not used to determine the 
surface brightness. The profile was corrected for back- 
ground and for exposure. As noted previously, the back- 
ground at large radii is quite uncertain due to the fore- 



ground emission by the North Polar Spur (§ 2), and we 
include a large systematic error in the background because 
of this. As discussed in § 2, this systematic error covered 
the range of NPS contributions from zero to all of the emis- 
sion at large radii on the S3 chip. Thus, the errors in the 
resulting surface brightness values and model fit parame- 
ters should include the full range of possible background 
contributions. 
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Fig. 7. — The profile of the unresolved, soft band (0.3-1 keV) 
emission as a function of the projected radius r. The short dashed 
line shows a de Vaucouleurs profile which fits the optical surface 
brightness of the galaxy. The long dashed curve is a beta model fit, 
while the solid curve is the best-fit assuming that the emission is 
the sum of emission which is proportional to the optical light plus 
a beta model. 

The observed profile is shown in Figure 7 as a function 
of the projected radius r. The observed soft X-ray profile 
is much broader than the optical profile for the galaxy. 
The short dashed line shows a de Vaucouleurs profile with 
an effective radius of r e ff = 72", which is the azimuthally- 
averaged value for the optical light. A fit to the profile 
with the X-ray surface brightness proportional to the op- 
tical surface brightness is completely unacceptable. We 
tried to fit the profile using the standard beta model 

, . 2"! -3/3+1/2 

Ix(a) = I . ( 3 ) 

where r c is the core radius. This provided a acceptable 
fit, which is shown as the long dashed curve in Figure 7. 
The value of (3 — 0.335 ±0.004 implies a rather fiat profile, 
although fairly flat profiles are found in many X-ray bright 
elliptical galaxies (e.g., Forman et al. 1985; Trinchieri, Fab- 
biano, & Canizares 1986). On the other hand, the core 
radius is required to be very small, r c — 2'.'9 ± 0'.'2. This 
implies that the X-ray surface brightness is nearly a power- 
law function of the radius, Ix oc r~ lm . 

The spectrum and colors of the unresolved emission in- 
dicate that a portion of this emission is due to unresolved 
stellar sources. The very small value of the core radius 
in the beta model fit also suggests that the emission near 
the center may have a significant contribution from un- 
resolved stellar sources with a cuspy de Vaucouleurs pro- 
file. Thus, we also tried fitting the surface brightness pro- 
file with the sum of a de Vaucouleurs profile and beta 
model, with the effective radius of the de Vaucouleurs pro- 
file fixed by the optical profile as before. This led to a 
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fit which was only marginally better than the single beta 
model (Ax 2 = —1.10 for one additional fitting param- 
eters), shown as the solid curve in Figure 7. The nor- 
malizations of the two components are roughly consistent 
with the decomposition of the unresolved emission into 
truly diffuse and discrete emission in Table 3 in the vari- 
ous regions based on X-ray colors. In this fit, the values 
of r c and /3 are both very large and poorly constrained, 
r c > 296" and /3 > 1.24. Obviously, these values are poor 
determined since the core radius is larger than the largest 
radius at which the surface brightness is determined. This 
suggests a rather flat surface brightness distribution for 
the truly diffuse emission out to a large radius. On the 
other hand, a constant surface brightness is not an ac- 
ceptable fit for the truly diffuse emission (A% 2 = +12.85 
for two fewer fitting parameters). The surface brightness 
due to diffuse emission is fit acceptably with a fairly flat 
power-law distribution with Ix oc r ~ - 55 . 

In conclusion, the surface brightness distribution of the 
diffuse emission is uncertain due to corrections for the un- 
resolved point sources and the background including the 
NPS. However, it has a radially declining surface bright- 
ness which is considerably flatter than that of the optical 
stellar emission from the galaxy. Note that this is true 
even if the diffuse emission is not corrected for the unre- 
solved point sources (Fig. 7), and that the very extended 
emission occurs at surface brightness levels which exceed 
those of the background and the NPS by approximately 



one order of magnitude. 

6. SPECTRAL ANALYSIS 

At the time when this analysis was done, there were con- 
siderable uncertainties in the soft X-ray spectral response 
of the S3 chip below 0.7 keV (e.g., Markevitch et al. 2000). 
We found significant residual deviations in most of the 
spectral fits for any reasonable spectrum at energies <0.7 
keV. Also, we have previous determined the spectrum of 
NGC 4697 with both the ROSAT PSPC and ASCA (ISB). 
We found that the Chandra spectrum at low energies was 
inconsistent with the ROSAT PSPC or the joint ROSAT 
and ASCA spectra. In addition to the spectra presented 
in ISB, we extracted the total ROSAT PSPC spectrum 
for the inner effective radius for comparison to the equiv- 
alent Chandra S3 spectrum for exactly the same region 
(§ 6.1). (We could not determine the ASCA spectrum of 
the same region as it is smaller than the spatial resolution 
of ASCA.) The Chandra and ROSAT spectra were com- 
pletely inconsistent at soft X-ray energies; they agree very 
well above 0.71 keV, and very poorly below this energy. 
Thus, we restrict all of our spectral analysis to the energy 
range 0.72-10.0 keV. In order to allow \ 2 statistics to be 
used, all of the spectra were grouped to at least 20 counts 
per spectral bin. 

Previous ROSAT and ASCA spectra of early- type galax- 
ies have indicated that they have at least two spectral 
components, a very hard component which may be due to 
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X-ray binaries and/or and AGN (Matsumoto et al. 1997; 
Allen et al. 2000), and a softer component. In X-ray lu- 
minous early-type galaxies, the soft component is domi- 
nant, and it is clearly due to diffuse gas at a temperature 
of ~1 keV (Forman et al. 1985; Canizares, Fabbiano, & 
Trinchieri 1987). In X-ray faint early- type galaxies, the 
soft component is much softer, and its origin is still uncer- 
tain (Fabbiano, Kim, & Trinchieri 1994; Pellegrini 1994; 
Kim et al. 1996; Irwin & Sarazin 1998a,b). The ASCA 
spectrum of the hard component has generally been fit 
by either a power-law (characterized by a photon spectral 
index T, where T > implies a photon spectrum which 
declines with energy; Allen et al. 2000) or by a thermal 
brcmsstrahlung spectrum (characterized by a hard com- 
ponent temperature Th\ Matsumoto et al. 1997). The soft 
component in X-ray bright galaxies is usually fit by the 
MEKAL model for the emission from a low density, opti- 
cally thin plasma (Irwin & Sarazin 1998a, b) This model 
is characterized by the temperature of the gas (T s ) and 
by the abundances of the heavy elements. Given the lim- 
ited statistics we have in our spectra, we will assume that 
the heavy element abundances have the solar ratios, and 
only allow the overall abundance of the heavy elements 
to vary. In X-ray faint galaxies, it is unclear what the 
appropriate soft emission model should be. If the soft 
emission is due to diffuse gas, then the MEKAL model 
would again be appropriate. If it is due to an optically 
thick stellar component (including the same LMXBs which 
produce the hard component), then it might be better 
represented as a blackbody, characterized by a temper- 
ature (T s again). Thus, we have used a spectral models 
which include both a hard (power-law or bremsstrahlung) 
and soft (MEKAL or blackbody) component. In ISB, the 
spectrum of NGC 4697 was fit with the sum of a hard 
bremsstrahlung and soft MEKAL model. 

Models were fit to the spectra using XSPEC. The re- 
sults are summarized in Table 4, where the errors are at 
the 90% confidence level. This Table gives the absorbing 
column Nh, hard spectral model, brcmsstrahlung temper- 
ature Th or power-law photon index T, unabsorbed flux 
of the hard component F\ (0.3-10 keV), soft component 
model, soft component temperature T s , abundance (rel- 
ative to solar) for the MEKAL component, unabsorbed 
flux of the soft component (0.3-10 keV), number of 
net counts in the spectrum, and \ 2 per degree of freedom 
(dof). The Galactic absorbing column towards NGC 4697 
is N H = 2.12 x 10 20 cm" 2 (Stark et al. 1992). 

6.1. X-ray Spectra of Resolved Sources 

The background for the resolved sources came from re- 
gions located around each source with three times the area 
of the regions for the source spectrum. We first determined 
the X-ray spectrum of all of the resolved sources within 
one effective radius. The resulting spectrum, which gave 
about 1600 counts in the 0.72-10 keV energy range after 
subtraction of background, is shown in Figure 8. 

We first fit the spectrum assuming the same spectra 
properties and properties as were used to fit the total 
ASCA plus ROS AT spectrum in ISB (row 1, Table 4). The 
normalizations of the hard (bremss) and soft (MEKAL) 
components were allowed to vary. This gave a good fit to 
the spectrum of the sources, but the flux of the soft com- 



ponent was very weak; the fit only gave an upper limit 
which was much lower than the hard flux. Thus, we tried 
removing the soft component completely in the fit (row 2). 
This gave a fit which was just as good (Ax 2 < 0.001) with 
one less free parameter. Thus, we conclude that, on av- 
erage, the resolved sources in the NGC 4697 have no soft 
component to their spectra, and that any soft component 
in the overall galaxy spectrum must come from a distinct, 
unresolved component, such as interstellar gas. 
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Fig. 8. — The data points in the upper panel are the cumulative 
X-ray spectrum of all of the resolved sources in the inner effective 
radius of NGC 4697. The solid histogram is the best-fit model spec- 
trum (row 3 in Table 4). The lower panel shows the residuals to the 
fit. 

Next, we allowed the temperature T/j of the hard compo- 
nent to vary (row 3). This resulted in a somewhat higher 
temperature than had been given by the total ASCA plus 
ROSAT spectrum, which did not resolve the spectra of 
the sources separately. The new fit was a considerable 
improvement, with A\ 2 = 8.7 for only one more free pa- 
rameter, which is significant at the 99.9% level according 
to the f-test. We will adopt this as our best-fit model for 
the resolved source population within the inner effective 
radius. This fit is shown with the data in Figure 8. Since 
this region combines good statistics with a high probabil- 
ity that the sources are associated with NGC 4697, we will 
also use this fit as our standard model for the spectrum of 
the LMXB population of NGC 4697. 

Next, we allowed the absorbing column to vary to see 
if there was evidence of any excess absorption beyond the 
Galactic column (row 4). We can only give an upper limit 
to the absorbing column, which is somewhat greater than 
the Galactic value. Presumably, this is partly due to re- 
moving the softest energies from the spectra due to the 
calibration uncertainties discussed above. Since allowing 
the absorbing column did not result in a greatly improved 
fit (A% 2 = 1.1 for one more free parameter), we will fix 
the absorption at the Galactic value. 

Finally, we considered different models for the hard and 
soft components. We represented the hard component by 
a power-law rather than by bremsstrahlung (row 5). This 
resulted in a slightly worse fit than with bremsstrahlung, 
so we will use the bremsstrahlung model for the hard com- 
ponent due to LMXBs. Then, we tried using a blackbody 
model for the soft component. This produced a better fit 
than the MEKAL model, but still no soft component was 
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required at the 90% confidence level. 

We also determined the cumulative X-ray spectrum of 
the resolved sources outside of one effective radius (rows 
7-9). As was true of the inner sources, these spectra are 
well-fit without any soft component in the spectrum. The 
cumulative spectra of the outer sources is considerably 
harder than that of the inner sources, and is fit best with 
an absorbing column which is much higher than Galactic 
(row 9). In row 7, we give the results of fitting the spec- 
trum of the outer sources using the same spectra parame- 
ters as fit the inner sources (row 3). This did not provide 
a very good fit. Rows 8 & 9 show the result of freeing 
the hard component temperature and the absorbing 
column. The inner source spectral model can be rejected 
with very high confidence (3>99%). We believe that this 
difference is mainly due to the greater influence of back- 
ground sources in the outer regions of the S3 chip. Many 
of the background sources have hardness ratios which in- 
dicate that they have very hard, strongly absorbed spectra 
(§ 4.4, Figure 4 in Paper I). For example, Src. 82 is the 
brightest X-ray source on the S3 chip, and provides about 
one quarter of the counts in the spectrum of outer sources. 
It has a very hard, strongly absorbed spectrum. 

We also tried to fit the spectra of the brightest and 
faintest sources separately (rows 10-13). Because the lu- 
minosity function shows a break at a characteristic lu- 
minosity Lb, we separated sources into brighter sources 
{L*x > Lb) and fainter sources (Lx < Lb). In order to 
include as many sources as possible to improve the statis- 
tics, but avoid too many background sources, we consid- 
ered all sources within an elliptical optical isophotc with a 
semimajor axis of 2 a c g . For the fainter sources, the best 
fit spectra were consistent with the Galactic column, so 
we fixed the absorption at this level. Unlike the case for 
the integrated spectra of all the sources within one effec- 
tive radius, the spectral fit was improved significantly by 
including a soft component. Either a MEKAL or black- 
body spectrum did almost exactly as well. However, the 
abundance in the MEKAL model was essentially uncon- 
strained, and there is no clear evidence for emission lines 
in the spectrum, so we adopt a blackbody model for the 
soft component. Rows 10 and 11 of Table 4 show the fits 
without and with the soft component. The flux of the soft 
component is small, but it improves the fit significantly. 
The spectra of the brighter sources did not require any 
soft component. On the other hand, they have a higher 
absorbing column and a lower bremsstrahlung tempera- 
ture. The net effect is that the fainter sources may have 
more very soft emission, while the brighter sources have 
more emission at 1-2 keV. 

6.2. X-ray Spectra of Unresolved Emission 

We also determined the X-ray spectra of the unresolved 
emission in NGC 4697. Because the unresolved emission 
has a low surface brightness, and the background is high 
and uncertain due to the superposed North Polar Spur 
(§ 2), we were only able to extract a useful spectrum for the 
central portion of the unresolved emission. Thus, we will 
only consider the unresolved X-ray spectrum from within 
the inner one effective radius; even in this region, the net 
counts in the spectrum after subtraction of background are 
only «700, which is marginal for spectral analysis. The 



spectrum was extracted from the inner one effective ra- 
dius region, excluding small regions around each of the 
sources. Because the excluded regions also contain un- 
resolved emission, the fluxes we derive will be somewhat 
smaller than the total fluxes for unresolved emission in the 
same region (e.g., as given in Table 3). The background 
for this spectrum and the spectra of the total emission 
(§ 6.3) were taken from a combination of blank sky and 
outer S3 chip spectra in order to remove the effects of the 
North Polar Spur, as discussed above (§ 2). We include 
a systematic error in the background, which means that 
the best-fit values of x 2 wm be somewhat smaller than 
expected for uncorrelated gaussian statistical errors. 
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Fig. 9. — The spectrum of unresolved emission from within the 
inner effective radius. The notation is the same as in Figure 8. 
The solid histogram is the ROSAT adapted spectral model for the 
unresolved emission (row 17 in Table 4). 

We started by fitting this spectrum with the spectrum 
derived from the ROSAT spectrum of the same region, 
which consisted of a bremsstrahlung hard component and 
a MEKAL soft component (row 14). Both the hard and 
soft component are required to fit the spectrum, but the 
hard component is fairly weak. The unresolved emission is 
much softer than the emission from the resolved sources. 
Allowing the absorbing column to vary did not improve the 
fit very significantly and resulted in only an upper limit 
on Nh (row 15). Thus, we will keep the column fixed at 
the Galactic value. 

If all of the spectral parameters are allowed to vary (row 
16), the fit is not significantly better, and the temperature 
of the hard bremsstrahlung component is nearly uncon- 
strained. Under the assumption that the hard component 
is from fainter discrete sources (presumably LMXBs) with 
a spectrum similar to that of the resolved sources, we fix 
the temperature of the hard component to the value (8.1 
keV) found for the resolved sources in the same region (row 

This fit is shown in row 17 and Figure 9. The observed 
spectrum shows the (mainly) Fe XVII (^0.72 keV emit- 
ted), Fe XVII (0.826 keV emitted), and the Ne X, Fe XVII, 
Fe XXI (~1.02 keV) line complexes, which indicates that 
the soft unresolved emission is due to diffuse gas. At lower 
energies where the spectral calibration is uncertain, the 
O VIII line at 0.654 keV is clearly seen. Some of the line 
ratios in the data appear to be inconsistent with the model; 
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for example, the Fe XVII (0.826 keV emitted) line complex 
is too weak in the model. This may indicate that there are 
calibration problems with the spectra even above 0.72 keV, 
or that the model is incorrect in some way. 

The fraction of the emission in the hard component is 
smaller than expected from the analysis of the X-ray col- 
ors (Table 3) or the spatial distribution of the unresolved 
emission (§ 5.3). This may be the result of uncertainties 
in the soft X-ray spectral response of the ACIS S3 de- 
tector, as noted above. Also, the regions used to exclude 
sources from the unresolved emission remove a higher frac- 
tion of the area near the center of the galaxy, and there 
is evidence that the unresolved emission gets softer as the 
radius increases. 

Total emission within 1 a_eff 
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Fig. 10. — The total (sources plus unresolved emission) spectrum 
from within the inner effective radius. The notation is the same as 
in Figure 8. The solid histogram is the ROSAT spectral model for 
the same region (row 18 in Table 4). 

6.3. Total X-ray Spectra 

We also determined the total (resolved and unresolved 
emission) X-ray spectra for various regions in NGC 4697. 
We first consider the total X-ray spectrum from within the 
inner one effective radius. Initially, we fit this spectrum 
using exactly the ROSAT spectrum (ISB) from the same 
region (row 18). This model included a bremsstrahlung 
hard component and a MEKAL soft component. As dis- 
cussed above, this gave a terrible fit if the energy range 
used was 0.3-10.0 keV, but gave a good fit if the soft X-ray 
channels were dropped (for the range 0.72-10.0 keV). The 
observed spectrum and the ROSAT model fit are shown 
in Figure 10. We then tried freeing the normalizations of 
the hard and soft components but leaving the shape of the 
two components fixed at the ROSAT form. This did not 
improve the fit, and the normalization were not changed 
significantly. We also tried allowing the absorbing column 
to vary (row 19); this did not improve the fit significantly, 
so we kept it fixed at the Galactic value. We then allowed 
all of the spectral parameters except Nh to vary (row 20); 
this did not provide a significantly improved fit given the 
five additional free parameters. Moreover, the abundance 
of the soft MEKAL emission was very poorly determined; 
since most of the soft emission in the spectrum is in lines, 
it is only the product of the MEKAL normalization times 
the abundance which is fixed by the spectrum. 

As was the case with the unresolved spectrum, the total 



spectrum shows lines due to (mainly) O VIII (0.654 keV 
emitted), Fe XVII (-0.72 keV emitted), Fe XVII (0.826 
keV emitted), and the Ne X, Fe XVII, Fe XXI (-1.02 keV) 
line complex. As was true of the unresolved spectrum, the 
Fe XVII (0.826 keV emitted) line complex is too weak in 
the model. Figure 10 also shows the Si Ka line at about 
1.8 keV in the data but not in the model. This is a back- 
ground feature produced by fluorescence in Chandra; its 
presence in the spectrum after subtraction of background 
indicates that the background subtraction is not perfect. 
The same feature is seen in the unresolved spectrum Fig- 
ure 9. 

At larger radii, the decreasing surface brightness of the 
emission, increasing contribution of the background, and 
the systematic uncertainty in the background make it more 
difficult to extract useful spectral information. We ex- 
tracted the spectrum for the region from 1 to 2 effective 
radii (row 21). This spectrum had only about 800 net 
counts after correction for the background. The tempera- 
ture of the hard component was not well determined, but 
was consistent with that found for the sources in the inner 
region (8.1 keV), so we adopt this value. The abundance 
in the soft component was not well-determined but con- 
sistent with the value from ROSAT and ASCA of 0.07, so 
we adopt this. The resulting spectral fit shows that the 
spectrum is more dominated by the soft emission than at 
smaller radii. There is no evidence for a significant tem- 
perature gradient, although the errors are large. 

7. DISCUSSION 

7.1. Nature of the Central Source 

The position of Src. 1 agrees with the optical position of 
the center of NGC 4697 to within the combined X-ray and 
optical errors. This suggests that this source is actually an 
AGN. The luminosity of this source, L x = 8 x 10 38 ergs 
s _1 , is higher than most LMXBs, but there are four sources 
which are brighter in the S3 field. Of these, the two bright- 
est sources in the field (Srcs. 72 & 82) are probably back- 
ground AGNs (Src. 72 definitely is), while the other two 
(Srcs. 36 & 40) are probably LMXBs in NGC 4697. On the 
other hand, the hardness ratios for the central source (H21 
= -0.22±°;°s, H31 = -0.34±°;°^) are very typical of those 
for the LMXBs (Table 2, Figure 4 in Paper I), which have 
average hardnesses ratios of (H21,H31) = (—0.14,-0.37) 
within the inner effective radius. Thus, based on the X- 
ray evidence alone, the central source might either be the 
active nucleus in NGC 4697, or one or more LMXBs seen 
in projection against the nucleus. The central source ap- 
pears to be slightly extended, which might indicate that it 
is the result of more than one point source. 

NGC 4697 is not a radio source at a fairly restrictive 
level (Birkinshaw & Davies 1985). We were unable to find 
any other clear evidence for nuclear activity in the litera- 
ture. This might argue that the central source is not due 
to an AGN. 

Stellar dynamical measurements indicate that 
NGC 4697 has a central black hole with a mass Mbh = 
1.6 x 10 8 Mq, assuming our adopted distance of 15.9 Mpc 
(Gebhardt et al. 2000). If we take the observed X-ray lu- 
minosity of the central source to be an upper limit to the 
luminosity of an associated active nucleus, then the AGN 
luminosity is < 4 x 10~ 8 of the Eddington luminosity of 
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the central black hole. 

7.2. The LMXB Population 
7.2.1. X-ray to Optical Luminosity Ratio of LMXB 's 

If one compares the discrete source (resolved and unre- 
solved LMXBs) X-ray luminosity from within one a c ff with 
the optical luminosity within the same region, one finds 
an X-ray-to-optical ratio for the LMXBs of Lx(LMXB, 
0.3-10 keV)/L B = 8.1 x 10 29 ergs s" 1 L B ^. For compar- 
ison, we found a value of (LMXB, 0.3-10 keY)/L B = 
7.2 x 10 29 erg s" 1 L% in the SO galaxy NGC 1553 (Blan- 
ton, Sarazin, & Irwin 2001), although the systematic er- 
rors are larger in that case because the observation only 
samples the top of the X-ray luminosity function. Chan- 
dra observations of other ellipticals appear to show larger 
variations in the X-ray to optical ratio (R. Mushotzky & 
L. Angclini, private communication). For the observed 
luminosity function in NGC 4697 (eq. [1]), most of the 
luminosity comes from sources near the break luminosity. 
In NGC 4697, there are only ^15 of these. Thus, part of 
the variation in the X-ray to optical ratios may be due to 
statistical fluctuations in the number of bright LMXBs, or 
to individual temporal variations in their X-ray emission. 

A recent XMM-Newton observation of the bulge of M31 
resolves ^90% of the emission into LMXBs within a 5' 
radius of the center (Shirey et al. 2001). Using our spec- 
tral model to convert to our passband, and determining 
the optical luminosity in this region from the surface pho- 
tometry of Walterbos & Kennicutt (1988), we find an 
X-ray to optical luminosity ratio of Lx(LMXB, 0.3-10 
keV)/Ls = 6.0 x 10 29 ergs s _1 L^^. Previously, Irwin 
& Sarazin (1998a) used ROSAT observations of the bulge 
of M31 (Primini, Forman, & Jones 1993; Supper et al. 
1997) to determine the X-ray luminosity. Assuming 90% 
of the emission is from LMXBs (based on the XMM result) 
and converting to our passband gives Lx (LMXB, 0.3-10 
keV)/Ls = 6.5 x 10 29 ergs s _1 Lg^, in very good agree- 
ment with the XMM result. However, the X-ray to optical 
ratio in the bulge of M31 seems to be about 35% smaller 
than that in NGC 4697. This may be due to the steep lu- 
minosity function in M31, which causes there to be no very 
bright LMXBs (§ 7.2.2). Sources with L x > L b contribute 
much of the emission in NGC 4697. 

Previous ROSAT observations of the bulge of the Sa 
galaxy NGC 1291 indicate an X-ray to optical ratio of 
L x (LMXB, 0.3-10 keV)/L B = 8.9 x 10 29 ergs s" 1 L B ^ 
(Irwin & Sarazin 1998a), if one converts from the ROSAT 
hard band of 0.52-2.02 keV using our best-fit source spec- 
trum. This agrees with the value for NGC 4697. We have 
confirmed this determination of the X-ray-to-optical ra- 
tio in NGC 1291 in a recent Chandra observation (Irwin, 
Sarazin, & Bregman 2001). However, about 30% of the 
source emission in NGC 1291 comes from the very bright 
central point source, which is probably a central AGN. If 
this central point source is removed to determine the X- 
ray luminosity due to LMXBs, then the X-ray-to-optical 
ratio in NGC 1291 is similar to that in the bulge of M31, 
and about 50% smaller than in NGC 4697. It appears that 
there are variations in the LMXB X-ray-to-optical ratios 
of early-type galaxies and spiral bulges. 



7.2.2. Luminosity Function 

The differential luminosity function of LMXBs in 
NGC 4697 is well- fit by a broken power law (eq. [1]) over 
the range of luminosities from 5 x 10 37 to 10 39 ergs s _1 . 
The break luminosity is L^ = 3.2+q' 9 x 10 38 ergs s _1 , which 
is similar to the Eddington luminosity for spherical accre- 
tion onto a 1.4 M Q neutron star (2 x 10 38 ergs s _1 for 
hydrogen accretion). Of course, in addition to the statisti- 
cal uncertainties, the value of Lb is affected by systematic 
uncertainties in the distance to NGC 4697 and in the con- 
version from counts to luminosity. This agreement sug- 
gests that the sources with luminosities above the break 
luminosity are accreting black holes, while those below the 
break are predominantly neutron stars. This would imply 
that NGC 4697 contains >15 luminous X-ray binary sys- 
tems containing black holes. If the more luminous of these 
systems (L x ~ 10 39 ergs s _1 ) are limited by the Eddington 
luminosity, they must contain fairly massive (M > 8M ) 
black holes. 

If the luminosity function were a broken power-law 
which curved up rather than down (an "elbow" rather than 
a "knee"), the break might be due to the simple superposi- 
tion of two power-law source distributions. Then, the ap- 
parent break in the luminosity function would just be the 
point where the two distributions crossed, and wouldn't 
necessarily have any physical meaning. However, such a 
superposition obviously cannot produce a broken power- 
law which curves downward (a knee), as is observed in 
NGC 4697. This requires that there be at least two distinct 
source populations above and below the break luminosity. 

We have found similar broken-power law luminosity 
functions in the SO galaxy NGC 1553 (Blanton et al. 2001), 
and the bulge of the early-type spiral NGC 1291 shows a 
sharp break right at the Eddington-break, with no bright 
non-nuclear sources above this break (Irwin et al. 2001). 
Apparently, similar broken-power-laws are found to fit the 
LMXB populations in X-ray bright ellipticals as well (R. 
Mushotzky & L. Angclini 2000, private communication). 
In all cases, the break luminosity agrees with the Edding- 
ton luminosity of a 1.4 Mq neutron star, to within the 
errors. Thus, this appears to be a universal feature of 
the LMXB populations in optically luminous, old stellar 
systems. 

In principle, this "Eddington-break" luminosity could 
be used as a distance estimator for galaxies. Margon 
& Ostriker (1973) originally proposed that X-ray sources 
would be limited by the Eddington luminosity and that 
this might be useful as a distance estimator. Also, the use 
of the peak luminosity of X-ray bursts as a distance indi- 
cator in our Galaxy had been suggested (e.g., van Paradijs 
1978), since Type I burst are also thought to be Eddington- 
limited. The Eddington luminosity depends on very sim- 
ple physics (electron scattering and gravity), and is de- 
termined solely by the mass of the accreting object and 
the mean mass per electron. In practice, the composi- 
tion dependence resolves itself into two interesting cases, 
hydrogen-dominated plasmas (the usual case) or heavy ele- 
ment dominated plasmas (helium or heavier), which might 
occur if the donor star had lost its outer hydrogen-rich en- 
velope. Although the physics of stellar core collapse is 
complex, the fact that most neutron stars in binaries have 
masses of approximately 1.4 M Q is presumably related to 
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the Chandrasehkar mass, which is also determined only by 
simple quantum mechanics, gravity, and slightly by com- 
position. Thus, the apparent universality of the break lu- 
minosity of LMXBs in early-type galaxies is understand- 
able in terms of very simple physics. 

Unfortunately, it isn't possible to determine the break 
luminosity in an individual galaxy with great statistical 
accuracy, because the number of bright LMXBs is very 
limited (about 15 in NGC 4697). For example, the sta- 
tistical errors in the flux associated with the break lumi- 
nosity in NGC 4697 imply that the distance could only be 
determined to about ±20% if the actual break luminos- 
ity were accurately known. Other distance estimators in 
current use give statistical errors of about 10% for early- 
type galaxies (e.g., Tonry et al. 2001). One situation where 
the Eddington break luminosity might be particularly use- 
ful would be in a group or cluster of galaxies, where one 
could determine the total luminosity function for all of the 
LMXBs in the early-type galaxies, and derive an average 
distance for the group or cluster. 

Prior to Chandra and XMM-Newton, the Population II 
system with the best-determined LMXB X-ray luminos- 
ity function was the bulge of M31 (Primini, Forman, & 
Jones 1993; Supper et al. 1997; Garcia ct al. 2000; Shirey 
et al. 2001). The recent XMM-Newton observation shows 
that the cumulative luminosity function of bulge sources 
in M31 is well-fit by a broken power-law, with an exponent 
of -0.47 ±0.03 for 36.2 < \ogL x < 37.4 and -1.79 ±0.26 
for 37.4 < logLx < 38.1 (Shirey et al. 2001). This applies 
to the 0.3-12 keV band, which is similar to the one we use 
(0.3-10 keV). These results are in good agreement with 
the luminosity function for the bulge of M31 previously 
determined by ROSAT (Primini, Forman, & Jones 1993; 
Supper ct al. 1997). The brightest source in the bulge of 
M31 has a luminosity of Lx ~ Lb, so we can't compare the 
upper end of the luminosity function in NGC 4697 with 
M31. On the other end, our source detection limit was 
\ogLx > 37.7, so we can't compare to the low luminosity 
end of the luminosity function of the bulge of M31. How- 
ever, we can compare luminosity functions in the rather 
narrow range 37.7 < \ogL x < 38.07 (0.3-10 keV). In this 
range, the exponent of the differential luminosity function 
(cq. [1]) in the bulge of M31 is a = 2.79 ± 0.26 (1-cr er- 
rors), while the exponent for NGC 4697 is ai = 1.29±qH 
(90% errors). Thus, the luminosity function in M31 de- 
clines much more steeply with luminosity. It seems that 
this rapid decline, combined with the much smaller optical 
luminosity of the bulge of M31, might account for the lack 
of brighter X-ray sources in M31. 

7.2.3. Globular Clusters and LMXBs 

We have identified 7 X-ray sources at radii of >1.'5 in 
NGC 4697 with candidate globular clusters. Because glob- 
ulars arc not resolved in ground based observations at the 
distance of NGC 4697, some of these candidate globulars 
might actually be background AGNs. Indeed, we have 
found that one of the candidate globular clusters (Hanes 
1977, cluster # 24, our Src. 72). is a background AGN 
at z = 0.696 (C. Mullis 2000, private communication). At 
present, there arc no lists of globular clusters in the cen- 
tral regions (<1.'5) of NGC 4697, where the diffuse optical 
brightness of the galaxy makes it difficult to identify glob- 



ular clusters in ground-based observations. Unfortunately, 
most of the X-ray sources (55 out of 90) are within 1.'5, and 
within this region very few of these sources are expected 
to be unrelated to NGC 4697. A Hubble Space Telescope 
study to detect globular cluster candidates in the central 
region of NGC 4697 would be very useful. 

Beyond 1'.5, roughly 20% of the X-ray sources are as- 
sociated with globular cluster candidates. On the other 
hand, the candidate globular clusters contain about 0.1% 
of the optical light of the galaxy in this region. Thus, the 
chance of an optical star being the donor in a LMXB is 
a factor of ^200 larger for stars in globular clusters than 
for field stars in NGC 4697. This indicates that globular 
clusters in this elliptical galaxy are a very hospitable en- 
vironment for active compact binaries, as is true in our 
own Galaxy (e.g., Hertz & Grindlay 1983; White, Nagase, 
& Parmar 1995). This is generally believed to result from 
stellar dynamical interactions in globular clusters, which 
can produce compact binary systems. 

One intriguing suggestion is that globular clusters might 
be the only location where very compact low mass binaries 
can evolve, and that all of the LMXBs in old stellar sys- 
tems may have been formed in globulars (White, Kulka- 
rni, & Sarazin 2001). Of course, not all of the LMXBs 
are presently found in globulars, either in the bulge of our 
Galaxy or in NGC 4697. The field LMXBs might have 
been ejected from globular clusters by kick velocities re- 
sulting from supcrnovae, by stellar dynamical processes, 
or by the dissolution of the globular due to tidal effects. 

Above, we have noted that there appear to be variations 
in the LMXB X-ray to optical luminosity ratio of early- 
type galaxies. If LMXBs all are produced in globulars, 
then the total luminosity of LMXBs might correlate better 
with the globular cluster population rather than the total 
optical luminosity. Since the specific frequency of glob- 
ular clusters (the number per optical luminosity) varies 
from galaxy to galaxy, this might explain the variation in 
the LMXB X-ray-to-optical ratio. If one compares the 
globular cluster population in NGC 4697 with the bulge 
of M31, NGC 4697 has 8^1 times as many globulars, and 
its specific globular cluster frequency Sn per unit optical 
luminosity is ^5 times higher (Harris 1991). The X-ray lu- 
minosity of LMXBs is about 5 times higher in NGC 4697, 
and the X-ray-to-optical ratio is about 1.35 times higher. 
Thus, the LMXB luminosity ratio is intermediate between 
the ratio of the optical luminosities and the ratio of the 
numbers of globular clusters. Obviously, there are too few 
galaxies with well-determined LMXB populations to al- 
low a really meaningful test at present, but there should 
be many more from Chandra and XMM-Newton within 
the next few years. One confusing issue is the fact that 
the luminosity function of LMXBs also appears to vary 
widely (§ 7.2.2), which requires something more complex 
than simple scaling of the entire population with optical 
luminosity or globular cluster population. 

In general, the specific frequency of globular clusters 
Sn increases from late to early-type spirals, from spirals 
to SOs, from SOs to ellipticals, and from normal giant el- 
lipticals to cDs (e.g., Harris 1991). If LMXBs all are born 
in globular clusters, then one might expect the specific 
frequency of LMXBs and their X-ray to optical ratio to 
increase in the same sequence. Already, Chandra observa- 
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tions may suggest an increase in the frequency of LMXBs 
in going from normal ellipticals to cDs (R. Mushotzky & 
L. Angclini 2000, private communication). 

7.3. Interstellar Gas Emission 

7.3.1. Nature of the Diffuse Emission 

The unresolved X-ray emission associated with 
NGC 4697 appears to consist of (at least) two compo- 
nents: unresolved LMXBs and truly diffuse emission. The 
X-ray colors of the unresolved emission, its X-ray spec- 
trum, and its spatial distribution are all consistent with 
this decomposition. The colors and spatial analysis indi- 
cate that «40% of the counts and ~61% of the luminosity 
is due to unresolved LMXBs. A similar result is found 
by extrapolating the luminosity function of the LMXBs 
to lower luminosities. The X-ray spectrum gives a smaller 
fraction of unresolved LMXBs, but this might be due 
to the uncertain soft response of the ACIS S3 or to the 
varying area of the galaxy removed in order to exclude 
sources. 

There are several arguments which indicate that the re- 
maining truly diffuse emission is due to hot interstellar 
gas. First, this emission has a much more extended radial 
distribution than the optical light in the galaxy (Fig. 7), 
which suggests that it is not from stellar sources. Sec- 
ond, the image of the diffuse emission (Fig. 2, and Fig. 2 
in Paper I) shows that it is somewhat irregular, with an 
L-shaped morphology. Such asymmetries are more easily 
produced in the distribution of gas (e.g., by ram pressure) 
than in the distribution of stellar sources. The diffuse 
emission is also much rounder than the highly elliptical 
distribution of optical light. Third, the diffuse emission 
is predominantly soft, with most of the emission below 1 
keV. Finally, the X-ray spectrum of the unresolved emis- 
sion shows low-ionization X-ray lines from ions such at 
O VIII and Fe XVII. This clearly implies that much of the 
soft emission is from an optically thin gas. 

The surface brightness distribution of this soft gaseous 
component is very broad and flat, and the total mass does 
not appear to converge within the S3 field of view. As 
a result, it is difficult to determine the total mass of the 
interstellar gas. This problem is exacerbated by the fore- 
ground emission from the North Polar Spur, which adds a 
(presumably) fairly flat emission component at a similar 
temperature. However, we have estimated the mass within 
a spherical region with a 75" radius, which corresponds to 
the average effective radius. The mass of interstellar gas 
is about 1.8 x 10 8 M within this region. This is much 
smaller than the masses of hot gas typically found in X- 
ray bright ellipticals (e.g., Forman, Jones, & Tucker 1985). 
However, the total mass is probably a factor of >3 larger 
given the flat surface brightness distribution. 

7.3.2. Origin of the Soft Component in X-ray Faint 
Ellipticals 

X-ray faint early-type galaxies have a very soft X-ray 
spectral component whose origin has been uncertain. One 
suggestion was that this emission came from the same 
LMXBs which produce the hard component (Irwin & 
Sarazin 1998a,b). Now, the LMXBs in NGC 4697 do pro- 
vide a significant part of the soft X-ray emission by simple 
virtue of the fact that they dominate the total emission. 



Certainly, some of the discrete X-ray sources in NGC 4697, 
the supersoft sources, do have a very soft X-ray spec- 
tral component. However, the cumulative spectrum of the 
sources does not show a significant soft X-ray spectral com- 
ponent, and is well-fit by a hard thermal brcmsstrahlung 
spectrum with a reasonably high temperature (kTh = 8.1 
keV; § 6.1). Similar results were found for the SO galaxy 
NGC 1553 (Blanton et al. 2001). While there are still 
considerable uncertainties in the soft X-ray spectral re- 
sponse of the ACIS detectors on Chandra, it appears very 
unlikely that the LMXBs can account for the very soft 
spectral component seen previously with ROS AT in either 
NGC 4697 or NGC 1553. By extension, it appears that 
the soft emission in X-ray faint ellipticals is not primarily 
due to LMXBs. 

The XMM-Newton spectra of the LMXBs in the bulge 
of M31 also shows that, taken together, they do not have a 
strong soft spectral component (Shirey et al. 2001). This 
confirms the ROS AT result of Borozdin & Pricdhorsky 
(2000), although Irwin & Bregman (1999) had reached the 
opposite conclusion. 

The soft spectral component might also be due to fainter 
stellar sources such as active M stars or RS CVn bina- 
ries (e.g., Pellegrini 1994). There were strong energetic 
and spectral arguments against any of these sources (Pel- 
legrini & Fabbiano 1994; Irwin & Sarazin 1998a). If the 
soft emission were due to faint stellar sources, it would be 
expected to follow the optical distribution in NGC 4697 
quite closely. In fact, the soft component is more spatially 
extended, rounder, and too irregular to be due to a large 
number of faint stellar sources (§ 5.3). 

The spatial distribution of the soft diffuse emission in 
both NGC 4697 and NGC 1553 implies that it is due to 
diffuse gas. Even more importantly, its X-ray spectrum 
is dominated by low ionization X-ray lines from ions like 
O VIII and Fe XVII. This clearly shows that the soft spec- 
tral component in these two early-type galaxies is due to 
diffuse interstellar gas. The recent XMM-Newton obser- 
vation of the central bulge of M31 also shows that there 
is diffuse emission with an emission-line spectrum (Shirey 
et al. 2001). Thus, it appears that the source of the soft 
component in early-type galaxies is cool (kT w 0.2 — 0.3 
keV) interstellar g suggested by Pellegrini & Fab- 

biano (1994). 

7.3.3. Physical State of the Hot Interstellar Gas 

The interstellar gas in NGC 4697 is rather cool (kT = 
0.29 keV) in comparison to X-ray bright ellipticals. A sim- 
ilar result was found for the SO galaxy NGC 1553 (Blanton 
et al. 2001), the bulge of the Sa galaxy NGC 1291 (Irwin 
et al. 2001), and for the diffuse gas in the bulge of M31 
(Shirey et al. 2001). NGC 4697 has a low velocity disper- 
sion for a bright elliptical, a = 165 km s _1 (Faber et. al. 
1989). On the other hand, this galaxy also has a very high 
rotation velocity for an elliptical (Peletier et al. 1990), so 
the velocity dispersion may not be representative of the 
full depth of the galactic potential. The observed gas tem- 
perature is consistent with the temperature vs. velocity 
dispersion (kT — a) correlation derived from ROS AT ob- 
servations of elliptical galaxies (Davis & White 1996), al- 
though this relationship was derived including a ROS AT 
single-temperature determination for NGC 4697 which is 
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higher than the Chandra temperature of the gas (presum- 
ably because LMXBs and gas could not be separated in 
the ROSAT spectrum). 

Why are the gas temperature and gas X-ray luminosity 
so low in NGC 4697 (and other X-ray faint ellipticals)? 
Pellegrini & Fabbiano (1994) argue that the low temper- 
ature is a result of a shallow gravitational potential, as 
indicated by the low velocity dispersion. As noted above, 
this is consistent with the observed temperature and a, al- 
though the strong rotation in this galaxy implies that the 
gravitational potential may be considerably deeper than 
expected from the velocity dispersion. 

What is the dynamical state of the interstellar gas in 
NGC 4697? X-ray bright ellipticals are believed to have 
global cooling inflows (e.g., Sarazin 1990). A global inflow 
can occur if the radiative cooling time of the gas is shorter 
than the age of the galaxy throughout the bulk of the op- 
tical image of the galaxy. Based on the best-fit model for 
the radial distribution of the soft X-ray emission, the cen- 
tral gas electron number density is about n e = 8.8 x 10 -3 
cm~ 3 . The central cooling time in the gas is about 4 x 10 8 
yr, which is much shorter than the likely age of the galaxy. 

The total X-ray luminosity of the gas may be too low 
to be consistent with a global cooling inflow involving the 
total stellar mass loss rate of the galaxy. The total bolo- 
metric luminosity of the gas is Lboi ~ 3.6 x 10 40 ergs s~ 4 , 
assuming the best-fit soft component model for the unre- 
solved emission (Table 4). Given its optical luminosity and 
assuming a normal elliptical galaxy stellar mass loss rate, 
the total rate of stellar mass loss is about 0.5 M Q yr -1 
(e.g., Sarazin 1990). The heat input to the gas which is 
just associated with the motions of the gas-losing stars is 
3cr 2 /2 per unit mass, which leads to a total heating rate of 
Eheat > 1-3 x 10 40 ergs s" 1 . If the gas involved in a cooling 
inflow, then infall in the galactic potential and adiabatic 
compression increase the heating rate by a factor of ~3. 
Type la supernova could also increase the heating rate by 
at least as large a factor. Thus, for a global cooling inflow 
one would expect a bolometric luminosity of > 8 x 10 40 
ergs s _1 . At best, this would be be marginally allowed by 
the observed luminosity. 

However, if the gas were involved in a global cooling 
flow, one would expect the gas X-ray surface brightness to 
be much more centrally peaked than is observed (Fig. 7). 
In global cooling flow models, the X-ray surface bright- 
ness is generally at least as centrally peaked as the optical 
surface brightness (e.g., Sarazin & White 1988; Sarazin & 
Ashe 1989). Thus, a global cooling flow seems unlikely. 
On the other hand, if the interstellar gas were in a global 
supersonic outflow or wind, the expected X-ray luminosity 
(< io 39 ergs s -1 ) would be much smaller than is observed 
(D'Ercolc ct al. 1989). 

Pellegrini & Fabbiano (1994) suggest that the gas in 
NGC 4697 is in a partial wind, with a cooling inflow at 
the center and an outflow in the outer regions. This could 
explain the low luminosity and gas mass in this galaxy. On 
the other hand, in these models the X-ray surface bright- 
ness of the gas is always more centrally condensed than the 
optical surface brightness (Pellegrini & Ciotti 1998), which 
is exactly the opposite of what we observe in NGC 4697 
(Fig. 7).' 

Among simple spherical hydrodynamical models of iso- 



lated galaxies, the only models which seem consistent with 
the low X-ray luminosity, low temperature, short cooling 
time, and broad distribution of the gas in NGC 4697 would 
seem to be subsonic inflation models prior to cooling in- 
flow (Ciotti et al. 1991; David, Forman, & Jones 1991). 
However, these models would probably require that much 
of the gas in the galaxy have been removed < 10 9 yr ago 
to explain the low mass of interstellar gas today, and also 
why this galaxy is at such a early hydrodynamical stage 
when it appears to have a very old stellar population. 

Pellegrini (1994) suggested that SO galaxies and other 
rapidly rotating early-type galaxies might be more suscep- 
tible to winds and other outflows as a result of rotation. 
She showed that there was an apparent anticorrclation be- 
tween the X-ray to optical ratio of early-type galaxies and 
their fractional rotational support. NGC 4697 might be 
an example of the effects of rapid rotation; this could help 
to explain the low mass of ISM gas. On the other hand, 
in this model one might expect the soft X-ray image to be 
highly elongated in the same direction as the optical im- 
age, which doesn't really seem to be the case (Figure 2). 

Alternatively, it may be that the gas in NGC 4697 has 
been strongly affected by ram pressure or other hydro- 
dynamical effects its environment (e.g., White & Sarazin 
1991). The extended morphology of the gas and its dis- 
torted morphology would suggest motion of the galaxy 
to the northwest relative to the ambient intergalactic 
medium. If the motion also has a significant component 
along the line of sight, ram pressure might explain both 
the extended distribution and the distorted morphology. 
Of course, ram pressure stripping would also explain the 
low X-ray luminosity and low gas mass. 

The best-fit abundance for the gas in NGC 4697 is rather 
low (0.06 solar), but is very poorly determined since much 
of the soft X-ray emission is due to lines. In any case, the 
abundance is consistent with the abundance-temperature 
correlation found by Davis & White (1996). 

7.3.4. Relation to Cooler Interstellar Medium 

NGC 4697 contains a small dust lane near its center 
(Goudfrooij et al. 1994; van Dokkum & Franx 1995), which 
is most prominent in images with the Hubble Space Tele- 
scope (e.g., Dejonghe et al. 1996). The dust lane is fairly 
regular, and appears to be part of a disk of material which 
is elongated along the optical major axis of the galaxy. 
There is no obvious corresponding feature in the X-ray 
image (Figure 1 and Paper I, Figure 2). 

NGC 4697 also contains an arc-shaped region of Ha 
emission, mainly extending to the northeast of the galaxy 
center (Goudfrooij et al. 1994). Again, there is nothing 
which obviously corresponds to this in the X-ray image 
(Figure 1 and Paper I, Figure 2). There does appear to be 
a bit more diffuse emission and more sources on this side of 
the galaxy, but this is not a very strong effect. The spatial 
decomposition of the unresolved emission suggests that it 
is dominated by unresolved LMXBs in the central parts of 
the galaxy. Thus, there may be too little hot gas to inter- 
act with the cooler ISM, and/or any spatial features due 
to interaction of diffuse X-ray gas with cooler ISM may be 
washed out by emission by unresolved LMXB emission. 

If the columns of cooler gas are significant, one might 
expect to see absorption features in the X-ray image cor- 
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responding to locations of cooler ISM. The lack of such 
features suggests that the columns of cooler gas are not 
very large (< 6 x 10 20 cm~ 2 ). 

8. CONCLUSIONS 

The origin of the X-ray emission in X-ray faint elliptical 
galaxies has been a mystery. Unlike X-ray bright early- 
type galaxies where the emission is dominated by interstel- 
lar gas with kT w 1 kcV, X-ray faint galaxies have both a 
very hard (~5 kcV) and very soft (—0.2 keV) spectral com- 
ponents. The hard component was thought to be due to 
LMXBs or an AGN. The origin of the very soft component 
was very uncertain; LMXBs (Irwin & Sarazin 1998a,b), 
cool interstellar gas (Pellegrini & Fabbiano 1994), or faint 
stellar sources, such as M stars and RS CVn binaries (e.g., 
Pellegrini 1994) have all been suggested. 

Our high spatial resolution Chandra observation of the 
X-ray faint elliptical galaxy NGC 4697 resolves most of the 
X-ray counts (61% within one effective radius) into point 
sources. A total of 90 individual sources are detected, of 
which —80 are low mass X-ray binaries (LMXBs) associ- 
ated with NGC 4697. Much of the emission is resolved 
even in the softest band. The dominance of LMXBs im- 
plies that this and other X-ray faint ellipticals have lost 
much of their interstellar gas. 

On the other hand, NGC 4697 does have a modest 
amount of X-ray emission by hot interstellar gas. Of the 
unresolved emission, it is likely that about half is from 
fainter LMXBs, while about 23% of the total count rate 
is due to interstellar gas with kT = 0.29 kcV. The X-ray 
spectrum of the soft diffuse emission is fit by a thermal 
model, and shows clear evidence for soft X-ray lines char- 
acteristic of low density, optically thin gas. The faint, 
diffuse soft X-ray emission in the Chandra image is very 
extended, indicating that the ISM gas distribution is much 
more extended and less centrally peaked than that of the 
optical light. The spatial distribution of the gas is rounder 
than that of the optical light, and has a somewhat irreg- 
ular, L-shaped form. This may indicate that the gas has 
been affected by ram pressure from motions through in- 
tergalactic gas. 

The X-ray luminosity and spatial distribution of the in- 
terstellar gas in NGC 4697 are inconsistent with a global 
spherical supersonic wind, partial wind, or global cooling 
flow. The gas may be undergoing subsonic inflation, or a 
rotationally induced outflux, or may be affected by ram 
pressure from intergalactic gas. There is no apparent re- 
lation between the X-ray emitting interstellar gas and the 
dust disk or Ha emission regions near the center of the 
galaxy. 

X-ray spectral analysis of the resolved and diffuse emis- 
sion indicates that the interstellar gas is the source of the 
very soft spectral component seen with ROSAT in the 
spectra of this and other X-ray faint early-type galaxies. 
Although the LMXBs produce a large amount of soft X-ray 
emission and some of them (the supersoft sources) do have 
a strong soft X-ray spectral component, taken together the 
spectrum of the LMXBs does not show a significant soft 
component. The cumulative LMXB spectrum is well- fit by 
thermal bremsstrahlung at kT = 8.1 kcV. The soft com- 
ponent in the spectrum is produced mainly by interstellar 
gas. 



NGC 4697 has an X-ray source located within 1" of the 
optical center with an X-ray luminosity of Lx — 8 x 10 38 
ergs s _1 (0.3-10 keV). This source may be due to an AGN 
and/or one or more LMXBs. Stellar dynamical measure- 
ments indicate that NGC 4697 has a massive central black 
hole; the central X-ray luminosity implies that this BH is 
emitting at < 4 x 10~ 8 of its Eddington luminosity. 

Three of the resolved sources are supersoft sources (e.g., 
Kahabka & van den Heuvel 1997). These are among the 
most distant and luminous supersoft sources known. 

Seven of the resolved sources in the outer parts of 
NGC 4697 (about 20%) are coincident with candidate 
globular clusters. This implies that globular clusters 
are a very favorable environment to form LMXBs; the 
same result is true of our Galaxy (e.g., Hertz & Grindlay 
1983). We discuss the possibility that all of the LMXBs in 
NGC 4697 (and other old stellar systems) were formed in 
globular clusters. An HST observation is needed to deter- 
mine the globular cluster population of the inner regions 
of NGC 4697 where most of the LMXBs are located. 

The X-ray-to-optical luminosity ratio for the LMXBs 
in NGC 4697 is L X (LMXB, 0.3-10 keV)/i B = 8.1 x 10 29 
ergs s _1 Lg\. This is about 35% higher than the value 
for the bulge of M31 and about 13% higher than the value 
we found for the SO galaxy NGC 1553. This suggests that 
there are variations in the X-ray-to-optical ratios early- 
type galaxies and spiral bulges. Since much of the X-ray 
luminosity in NGC 4697 comes from a small number (—15) 
of bright sources, part of this variation may be due to 
statistical fluctuations in the number of bright sources or 
temporal variations in their individual emission. If most 
LMXBs are formed in globular clusters, then their num- 
ber may correlate more directly with the population of 
globulars. On the other hand, in the range they have 
in common, the luminosity function of the bulge of M31 
declines more rapidly with increasing luminosity than in 
NGC 4697. This suggests that the populations of LMXBs 
cannot simply be scaled from galaxy to galaxy. 

The X-ray luminosities (0.3-10 keV) of the resolved 
LMXBs range from -5 x 10 37 to -2.5 x 10 39 ergs s" 1 . The 
luminosity function has a "knee" at 3.2 x 10 38 ergs s _1 , 
which is approximately the Eddington luminosity of a 1.4 
Mq neutron star (NS). Based on other observations with 
Chandra, this knee appears to be a characteristic feature of 
the luminosity functions of LMXBs in early-type galaxies. 
This knee might provide a standard candle which could 
be used to determine distances to galaxies. This knee may 
separate accreting NS and black hole (BH) binaries. If 
the brightest sources in NGC 4697 are Eddington limited, 
they must contain fairly massive BHs. 

These and other Chandra observations of luminous el- 
liptical galaxies may represent the first direct detections of 
neutron stars and stellar mass black holes in these galax- 
ies. Our detection of a large population of binaries with 
NSs and massive BHs provides perhaps the most direct ev- 
idence that this elliptical galaxy (or its progenitors) once 
contained a large number of massive main sequence stars. 
The population of LMXBs provides a tool to study the 
high mass end of the initial mass function of early-type 
galaxies, long after the massive main sequence stars have 
died. 
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